(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 14.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    213871,       5226]
NotebookOptionsPosition[    202178,       5005]
NotebookOutlinePosition[    202617,       5022]
CellTagsIndexPosition[    202574,       5019]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Model + Two-sided Loyal Users", "Chapter",
 CellChangeTimes->{{3.9625997412273245`*^9, 3.962599764101927*^9}, {
  3.971341551694647*^9, 
  3.9713415569781494`*^9}},ExpressionUUID->"2f76dbf6-4856-fb4a-b09f-\
e7ae0b559ca0"],

Cell[CellGroupData[{

Cell["Utility Functions", "Section",
 CellChangeTimes->{{3.962592639745962*^9, 
  3.9625926487241325`*^9}},ExpressionUUID->"075a57bb-5672-684f-b1a0-\
fcf1d8427ebf"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"u", "[", 
    RowBox[{
    "\[Sigma]_", ",", "\[Alpha]_", ",", "n_", ",", "p_", ",", "t_", ",", 
     "m_"}], "]"}], ":=", 
   RowBox[{"\[Sigma]", "+", 
    RowBox[{"\[Alpha]", "*", "n"}], "-", "p", "-", 
    RowBox[{"t", "*", "m"}]}]}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "\[Alpha]", " ", "is", " ", "\[Alpha]", " ", "or", " ", "\[Beta]"}], ",", 
    " ", 
    RowBox[{"n", " ", "is", " ", "quantities"}], ",", 
    RowBox[{"t", " ", "is", " ", "z", " ", "or", " ", "q"}]}], "*)"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.961794882527809*^9, 3.961794891489977*^9}, {
  3.961794986475697*^9, 3.9617951599176235`*^9}, {3.9617957652547646`*^9, 
  3.9617958228781624`*^9}, {3.962600117341297*^9, 3.962600163803043*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"670565b0-2745-7940-9f25-1ce8412da6af"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"u", "[", 
   RowBox[{"\[Sigma]x", ",", "\[Alpha]", ",", 
    RowBox[{"yA", "+", "LyA"}], ",", "pA", ",", "z", ",", "x"}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"u", 
    RowBox[{"(", 
     RowBox[{"x", ";", "A"}], ")"}]}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.961795169149374*^9, 3.961795171046116*^9}, {
  3.961795835469267*^9, 3.961795846230795*^9}, {3.962044354229542*^9, 
  3.9620443687135353`*^9}, {3.962044407542078*^9, 3.962044412323456*^9}, {
  3.962044496089268*^9, 3.962044496318153*^9}, {3.962044913496469*^9, 
  3.9620449275263844`*^9}, {3.962048146390795*^9, 3.962048151120805*^9}, {
  3.962051903154426*^9, 3.9620519045033784`*^9}},
 CellLabel->"In[2]:=",ExpressionUUID->"3109ff47-9610-f24b-b339-93226a554405"],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "pA"}], "-", 
  RowBox[{"x", " ", "z"}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"LyA", "+", "yA"}], ")"}], " ", "\[Alpha]"}], "+", 
  "\[Sigma]x"}]], "Output",
 CellChangeTimes->{
  3.96179517242914*^9, 3.961795848885988*^9, 3.9619693254139347`*^9, 
   3.962043442996704*^9, 3.9620444162913933`*^9, 3.962044499241476*^9, {
   3.962044892592314*^9, 3.96204492952683*^9}, 3.96204835450165*^9, 
   3.962051927219168*^9, 3.9621353193336735`*^9, 3.9621375094244556`*^9, 
   3.962592758207905*^9, 3.962598294834238*^9, 3.962598715591514*^9, 
   3.962598976798109*^9, 3.96259932719295*^9, 3.962599951525444*^9, 
   3.9713410632251015`*^9},
 CellLabel->"Out[2]=",ExpressionUUID->"fd18633c-2e83-084a-9f4d-a74f97e793b4"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"u", "[", 
   RowBox[{"\[Sigma]x", ",", "\[Alpha]", ",", 
    RowBox[{"Y", "-", "yB", "+", "LyB"}], ",", "pB", ",", "z", ",", 
    RowBox[{"1", "-", "x"}]}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"u", 
    RowBox[{"(", 
     RowBox[{"x", ";", "B"}], ")"}]}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.961795127356207*^9, 3.9617951300439777`*^9}, {
   3.961795180283188*^9, 3.9617952446505814`*^9}, 3.961795394517208*^9, {
   3.96179585270656*^9, 3.961795858210554*^9}, {3.962044422333343*^9, 
   3.962044430102932*^9}, {3.962044502578022*^9, 3.962044502832964*^9}, {
   3.962044838705118*^9, 3.962044847225445*^9}, {3.9620449384214554`*^9, 
   3.962044943318269*^9}, {3.9620481615964823`*^9, 3.962048164778518*^9}, {
   3.962048359898278*^9, 3.96204836518503*^9}, {3.9620519076724386`*^9, 
   3.9620519088784366`*^9}},
 CellLabel->"In[3]:=",ExpressionUUID->"ecf37400-b462-1648-83b7-dc98d91d3794"],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "pB"}], "-", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"1", "-", "x"}], ")"}], " ", "z"}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"LyB", "+", "Y", "-", "yB"}], ")"}], " ", "\[Alpha]"}], "+", 
  "\[Sigma]x"}]], "Output",
 CellChangeTimes->{3.961795247815909*^9, 3.961795862035525*^9, 
  3.961969328654583*^9, 3.962043445704848*^9, 3.962044433976755*^9, 
  3.9620445051284504`*^9, 3.9620448562034264`*^9, 3.962044894835424*^9, 
  3.9620449453609047`*^9, 3.962048368352512*^9, 3.96205192734069*^9, 
  3.9621353193694553`*^9, 3.962137509480511*^9, 3.962592758237924*^9, 
  3.9625982948647423`*^9, 3.96259871562562*^9, 3.9625989768545437`*^9, 
  3.9625993272358418`*^9, 3.9625999515623856`*^9, 3.9713410632764893`*^9},
 CellLabel->"Out[3]=",ExpressionUUID->"152e8bff-350e-bb4b-ab1c-2c3dfae3c5ef"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"u", "[", 
   RowBox[{"\[Sigma]y", ",", "\[Beta]", ",", 
    RowBox[{"xA", "+", "LxA"}], ",", "rA", ",", "q", ",", "y"}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"u", 
    RowBox[{"(", 
     RowBox[{"y", ";", "A"}], ")"}]}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.9617953470045795`*^9, 3.96179548090883*^9}, {
  3.961795512361624*^9, 3.961795565151766*^9}, {3.961795876481941*^9, 
  3.9617959311183243`*^9}, {3.961796007903618*^9, 3.961796031901535*^9}, {
  3.962044437479868*^9, 3.962044444875263*^9}, {3.962044507670639*^9, 
  3.9620445128377876`*^9}, {3.9620449495767384`*^9, 3.9620449598760376`*^9}, {
  3.962045010156624*^9, 3.962045015069948*^9}, {3.962048177138963*^9, 
  3.962048184162804*^9}, {3.962048375964861*^9, 3.962048378918632*^9}, {
  3.962051911591855*^9, 3.962051912787611*^9}},
 CellLabel->"In[4]:=",ExpressionUUID->"e73ed0d8-3fe4-ca40-82c3-d5ac7f8f9886"],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "rA"}], "-", 
  RowBox[{"q", " ", "y"}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"LxA", "+", "xA"}], ")"}], " ", "\[Beta]"}], "+", 
  "\[Sigma]y"}]], "Output",
 CellChangeTimes->{
  3.961796034706957*^9, 3.961969330983097*^9, 3.9620434484812584`*^9, 
   3.962044447050604*^9, 3.962044514945812*^9, {3.962044957642666*^9, 
   3.962044961370577*^9}, 3.962045016360571*^9, 3.962048381739895*^9, 
   3.962051927421889*^9, 3.962135319414879*^9, 3.9621375095046654`*^9, 
   3.962592758245479*^9, 3.962598294875105*^9, 3.9625987156661415`*^9, 
   3.9625989769011917`*^9, 3.9625993272721214`*^9, 3.96259995159639*^9, 
   3.971341063323969*^9},
 CellLabel->"Out[4]=",ExpressionUUID->"6ec3e942-74f5-db4a-99f1-8f8c9cecb68b"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"u", "[", 
   RowBox[{"\[Sigma]y", ",", "\[Beta]", ",", 
    RowBox[{"1", "-", "xB", "+", "LxB"}], ",", "rB", ",", "q", ",", 
    RowBox[{"Y", "-", "y"}]}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"u", 
    RowBox[{"(", 
     RowBox[{"y", ";", "B"}], ")"}]}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.961796043448314*^9, 3.9617960819692535`*^9}, {
  3.9620445182262383`*^9, 3.962044528101897*^9}, {3.962045023289093*^9, 
  3.962045029649912*^9}, {3.962048205121745*^9, 3.9620482083529835`*^9}, {
  3.962048387193147*^9, 3.962048388959099*^9}, {3.962051914947723*^9, 
  3.9620519154373813`*^9}},
 CellLabel->"In[5]:=",ExpressionUUID->"99bdc6d7-0dc0-5946-bfdd-8369c67bc491"],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "rB"}], "-", 
  RowBox[{"q", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "y"}], "+", "Y"}], ")"}]}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"1", "+", "LxB", "-", "xB"}], ")"}], " ", "\[Beta]"}], "+", 
  "\[Sigma]y"}]], "Output",
 CellChangeTimes->{3.961796084952175*^9, 3.961969333025295*^9, 
  3.962043450734501*^9, 3.962044530141474*^9, 3.962045031180122*^9, 
  3.9620483923274565`*^9, 3.962051927497985*^9, 3.962135319460934*^9, 
  3.962137509537443*^9, 3.9625927582707214`*^9, 3.962598294900778*^9, 
  3.962598715707821*^9, 3.962598976937456*^9, 3.962599327282217*^9, 
  3.962599951627104*^9, 3.971341063336273*^9},
 CellLabel->"Out[5]=",ExpressionUUID->"1a5004df-7251-8a4e-b222-4cf13462a382"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"u", "[", 
   RowBox[{
    RowBox[{"2", "\[Sigma]x"}], ",", "\[Alpha]", ",", 
    RowBox[{"Y", "+", "LyA", "+", "LyB"}], ",", 
    RowBox[{"pA", "+", "pB"}], ",", "z", ",", "1"}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"u", 
    RowBox[{"(", 
     RowBox[{"x", ";", 
      RowBox[{
       RowBox[{"A", "&"}], "B"}]}], ")"}]}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.961796093827093*^9, 3.96179613146328*^9}, {
  3.9617961993631077`*^9, 3.961796234055187*^9}, {3.962048243071376*^9, 
  3.962048252188691*^9}, {3.962048398442818*^9, 3.962048410363828*^9}, {
  3.962051947284481*^9, 3.9620519521641426`*^9}},
 CellLabel->"In[6]:=",ExpressionUUID->"5d32904f-bbb5-8b4f-bb36-c2fa067cb57b"],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "pA"}], "-", "pB", "-", "z", "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"LyA", "+", "LyB", "+", "Y"}], ")"}], " ", "\[Alpha]"}], "+", 
  RowBox[{"2", " ", "\[Sigma]x"}]}]], "Output",
 CellChangeTimes->{
  3.96179623635656*^9, 3.961969335747751*^9, 3.962043452869201*^9, 
   3.962044593691204*^9, 3.962048415160576*^9, {3.96205193655367*^9, 
   3.962051965403107*^9}, 3.962135319501024*^9, 3.962137509584439*^9, 
   3.9625927582987404`*^9, 3.962598294934778*^9, 3.962598715716789*^9, 
   3.962598976948973*^9, 3.962599327311769*^9, 3.962599951665081*^9, 
   3.971341063387272*^9},
 CellLabel->"Out[6]=",ExpressionUUID->"4c00b93b-502b-9f4b-a4ec-24319eb8e863"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"u", "[", 
   RowBox[{
    RowBox[{"2", "\[Sigma]y"}], ",", "\[Beta]", ",", 
    RowBox[{"1", "+", "LxA", "+", "LxB"}], ",", 
    RowBox[{"rA", "+", "rB"}], ",", "q", ",", "Y"}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"u", 
    RowBox[{"(", 
     RowBox[{"y", ";", 
      RowBox[{
       RowBox[{"A", "&"}], "B"}]}], ")"}]}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.961796251764866*^9, 3.961796267159481*^9}, {
  3.9617963049741325`*^9, 3.9617963666747875`*^9}, {3.962048430820587*^9, 
  3.962048435478168*^9}, {3.962051957254796*^9, 3.962051961048887*^9}},
 CellLabel->"In[7]:=",ExpressionUUID->"92b8c5ca-2962-cc49-b1d3-b66eb82809d2"],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "rA"}], "-", "rB", "-", 
  RowBox[{"q", " ", "Y"}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"1", "+", "LxA", "+", "LxB"}], ")"}], " ", "\[Beta]"}], "+", 
  RowBox[{"2", " ", "\[Sigma]y"}]}]], "Output",
 CellChangeTimes->{
  3.961796368967282*^9, 3.961969337822438*^9, 3.9620434554829597`*^9, 
   3.962044595713356*^9, 3.962048438915116*^9, {3.962051936724155*^9, 
   3.962051965537991*^9}, 3.9621353195415287`*^9, 3.962137509585785*^9, 
   3.962592758326071*^9, 3.962598294961796*^9, 3.9625987157490597`*^9, 
   3.962598976997547*^9, 3.962599327338663*^9, 3.962599951700676*^9, 
   3.9713410634386044`*^9},
 CellLabel->"Out[7]=",ExpressionUUID->"314ef302-c208-0942-96a4-07c67cd0bb57"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Indifferent Consumers", "Section",
 CellChangeTimes->{{3.9625926557625217`*^9, 
  3.9625926638043785`*^9}},ExpressionUUID->"3aa77a04-0892-7645-acba-\
1d32f23a05b0"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"xAhat", "[", 
   RowBox[{
   "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_", 
    ",", "pB_", ",", "z_", ",", "yBhat_", ",", "yAhat_", ",", "q_", ",", 
    "\[Beta]_", ",", "rA_", ",", "rB_", ",", "LyA_", ",", "LyB_"}], "]"}], 
  "=", 
  RowBox[{"xAhat", "/.", 
   RowBox[{
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"u", "[", 
         RowBox[{
          RowBox[{"2", "\[Sigma]x"}], ",", "\[Alpha]", ",", 
          RowBox[{"Y", "+", "LyA", "+", "LyB"}], ",", 
          RowBox[{"pA", "+", "pB"}], ",", "z", ",", "1"}], "]"}], "-", 
        RowBox[{"u", "[", 
         RowBox[{"\[Sigma]x", ",", "\[Alpha]", ",", 
          RowBox[{"Y", "-", "yBhat", "+", "LyB"}], ",", "pB", ",", "z", ",", 
          RowBox[{"1", "-", "xAhat"}]}], "]"}]}], "==", "0"}], ",", "xAhat"}],
      "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.961796388232682*^9, 3.961796408830221*^9}, {
   3.961796439708488*^9, 3.961796541795967*^9}, {3.9617967130632973`*^9, 
   3.961796749105871*^9}, {3.961796781846386*^9, 3.9617968994140625`*^9}, {
   3.961796931069947*^9, 3.9617969510105495`*^9}, {3.961797692635395*^9, 
   3.961797702564446*^9}, {3.9620440374700584`*^9, 3.96204405612393*^9}, {
   3.962044165378381*^9, 3.962044166575917*^9}, {3.962044247841337*^9, 
   3.962044248601233*^9}, {3.9620445599276905`*^9, 3.962044569965278*^9}, {
   3.962044602615799*^9, 3.962044614624906*^9}, 3.962045061183962*^9, {
   3.962048464787542*^9, 3.962048484557157*^9}, {3.96204851519029*^9, 
   3.962048558994462*^9}, {3.9620519867224007`*^9, 3.96205200663429*^9}},
 CellLabel->"In[8]:=",ExpressionUUID->"4aeeeca2-80de-ef43-b694-c81a806f132d"],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"-", "pA"}], "+", 
   RowBox[{"LyA", " ", "\[Alpha]"}], "+", 
   RowBox[{"yBhat", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], "z"]], "Output",
 CellChangeTimes->{
  3.961796967098982*^9, 3.961797706075032*^9, 3.961969341140175*^9, 
   3.9620434589515114`*^9, {3.962044585416033*^9, 3.9620445981886463`*^9}, 
   3.962044637141823*^9, 3.9620450799763527`*^9, 3.96204849408097*^9, 
   3.962048526889663*^9, {3.962051972414913*^9, 3.962052010987238*^9}, 
   3.962135319582306*^9, 3.962137509629566*^9, 3.962592758355032*^9, 
   3.962598294998228*^9, 3.962598715789465*^9, 3.962598983186827*^9, 
   3.962599332048958*^9, 3.962599951737902*^9, 3.9713410635292377`*^9},
 CellLabel->"Out[8]=",ExpressionUUID->"852cecc7-adc5-4340-83d6-2e81ee5c8ad2"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"xBhat", "[", 
   RowBox[{
   "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_", 
    ",", "pB_", ",", "z_", ",", "yBhat_", ",", "yAhat_", ",", "q_", ",", 
    "\[Beta]_", ",", "rA_", ",", "rB_", ",", "LyA_", ",", "LyB_"}], "]"}], 
  "=", 
  RowBox[{"xBhat", "/.", 
   RowBox[{
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"u", "[", 
         RowBox[{
          RowBox[{"2", "\[Sigma]x"}], ",", "\[Alpha]", ",", 
          RowBox[{"Y", "+", "LyA", "+", "LyB"}], ",", 
          RowBox[{"pA", "+", "pB"}], ",", "z", ",", "1"}], "]"}], "-", 
        RowBox[{"u", "[", 
         RowBox[{"\[Sigma]x", ",", "\[Alpha]", ",", 
          RowBox[{"yAhat", "+", "LyA"}], ",", "pA", ",", "z", ",", "xBhat"}], 
         "]"}]}], "==", "0"}], ",", "xBhat"}], "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.961797093560528*^9, 3.961797193183826*^9}, {
   3.961797579475991*^9, 3.9617975819408245`*^9}, 3.961797712897291*^9, {
   3.962045123710243*^9, 3.962045149201183*^9}, {3.962048576580638*^9, 
   3.9620486084236937`*^9}, {3.962052028684332*^9, 3.962052043755837*^9}},
 CellLabel->"In[9]:=",ExpressionUUID->"d716037c-e532-c942-844f-e2904edaa450"],

Cell[BoxData[
 FractionBox[
  RowBox[{"pB", "+", "z", "-", 
   RowBox[{"LyB", " ", "\[Alpha]"}], "-", 
   RowBox[{"Y", " ", "\[Alpha]"}], "+", 
   RowBox[{"yAhat", " ", "\[Alpha]"}], "-", "\[Sigma]x"}], "z"]], "Output",
 CellChangeTimes->{3.961797201411122*^9, 3.9617975854305096`*^9, 
  3.961797715439348*^9, 3.961969343477028*^9, 3.962043461823601*^9, 
  3.962045157385669*^9, 3.962048621924969*^9, 3.962052047628209*^9, 
  3.962135319618299*^9, 3.962137509685847*^9, 3.962592758385144*^9, 
  3.962598295028059*^9, 3.962598715826674*^9, 3.9625989834757977`*^9, 
  3.9625993323094425`*^9, 3.962599951770954*^9, 3.971341063640026*^9},
 CellLabel->"Out[9]=",ExpressionUUID->"5a157eb2-1602-d34d-a08c-ad7534086220"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"yAhat", "[", 
   RowBox[{
   "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_", 
    ",", "pB_", ",", "z_", ",", "xAhat_", ",", "xBhat_", ",", "q_", ",", 
    "\[Beta]_", ",", "rA_", ",", "rB_", ",", "LxA_", ",", "LxB_"}], "]"}], 
  "=", 
  RowBox[{"yAhat", "/.", 
   RowBox[{
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"u", "[", 
         RowBox[{
          RowBox[{"2", "\[Sigma]y"}], ",", "\[Beta]", ",", 
          RowBox[{"1", "+", "LxA", "+", "LxB"}], ",", 
          RowBox[{"rA", "+", "rB"}], ",", "q", ",", "Y"}], "]"}], "-", 
        RowBox[{"u", "[", 
         RowBox[{"\[Sigma]y", ",", "\[Beta]", ",", 
          RowBox[{"1", "-", "xBhat", "+", "LxB"}], ",", "rB", ",", "q", ",", 
          RowBox[{"Y", "-", "yAhat"}]}], "]"}]}], "==", "0"}], ",", "yAhat"}],
      "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.96179723659923*^9, 3.961797237867476*^9}, {
   3.961797311389097*^9, 3.9617973331604958`*^9}, {3.9617973965905857`*^9, 
   3.961797464839428*^9}, {3.961797571040652*^9, 3.9617976144991627`*^9}, 
   3.9617977183691*^9, 3.961797996115477*^9, {3.962045169645757*^9, 
   3.9620452120944424`*^9}, {3.962048644998523*^9, 3.962048686178957*^9}, {
   3.9620520646809597`*^9, 3.962052074611221*^9}},
 CellLabel->"In[10]:=",ExpressionUUID->"0f77b9f7-6d6a-d24f-9e59-8212aea8fc7b"],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"-", "rA"}], "+", 
   RowBox[{"LxA", " ", "\[Beta]"}], "+", 
   RowBox[{"xBhat", " ", "\[Beta]"}], "+", "\[Sigma]y"}], "q"]], "Output",
 CellChangeTimes->{
  3.961797669236237*^9, 3.961797720664686*^9, 3.9617979977254295`*^9, 
   3.9619693461901207`*^9, 3.962043464462049*^9, {3.962045200707926*^9, 
   3.962045214482681*^9}, 3.962048691048114*^9, 3.962052077714184*^9, 
   3.962135319664024*^9, 3.9621375097252846`*^9, 3.962592758411989*^9, 
   3.9625982950633965`*^9, 3.962598715860338*^9, 3.9625989835112267`*^9, 
   3.9625993323462906`*^9, 3.962599951807738*^9, 3.9713410637353*^9},
 CellLabel->"Out[10]=",ExpressionUUID->"e704a5ef-27b6-944c-afce-32ebcbbde64e"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"yBhat", "[", 
   RowBox[{
   "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_", 
    ",", "pB_", ",", "z_", ",", "xAhat_", ",", "xBhat_", ",", "q_", ",", 
    "\[Beta]_", ",", "rA_", ",", "rB_", ",", "LxA_", ",", "LxB_"}], "]"}], 
  "=", 
  RowBox[{"yBhat", "/.", 
   RowBox[{
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"u", "[", 
         RowBox[{
          RowBox[{"2", "\[Sigma]y"}], ",", "\[Beta]", ",", 
          RowBox[{"1", "+", "LxA", "+", "LxB"}], ",", 
          RowBox[{"rA", "+", "rB"}], ",", "q", ",", "Y"}], "]"}], "-", 
        RowBox[{"u", "[", 
         RowBox[{"\[Sigma]y", ",", "\[Beta]", ",", 
          RowBox[{"xAhat", "+", "LxA"}], ",", "rA", ",", "q", ",", "yBhat"}], 
         "]"}]}], "==", "0"}], ",", "yBhat"}], "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.961798009195709*^9, 3.961798068187819*^9}, {
   3.962045233908418*^9, 3.962045247934807*^9}, {3.962048700878498*^9, 
   3.962048709670955*^9}, 3.962048741575404*^9, 3.96204891075371*^9, {
   3.96205210758554*^9, 3.962052119448369*^9}},
 CellLabel->"In[11]:=",ExpressionUUID->"fac2a9da-1722-e446-b9f4-1085badaa6db"],

Cell[BoxData[
 FractionBox[
  RowBox[{"rB", "+", 
   RowBox[{"q", " ", "Y"}], "-", "\[Beta]", "-", 
   RowBox[{"LxB", " ", "\[Beta]"}], "+", 
   RowBox[{"xAhat", " ", "\[Beta]"}], "-", "\[Sigma]y"}], "q"]], "Output",
 CellChangeTimes->{3.9620460692827797`*^9, 3.962048754864609*^9, 
  3.962048914580781*^9, 3.962052122800015*^9, 3.962135319719225*^9, 
  3.96213750977297*^9, 3.962592758438034*^9, 3.962598295095859*^9, 
  3.962598715900749*^9, 3.962598983548525*^9, 3.962599332382572*^9, 
  3.9625999518417416`*^9, 3.9713410638309155`*^9},
 CellLabel->"Out[11]=",ExpressionUUID->"081f3fdd-ed1f-a646-a89b-4078e572e193"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{"xAhatSol", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_",
       ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", 
      "rB_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], ",", 
    RowBox[{"xBhatSol", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_",
       ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", 
      "rB_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], ",", 
    RowBox[{"yAhatSol", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_",
       ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", 
      "rB_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], ",", 
    RowBox[{"yBhatSol", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_",
       ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", 
      "rB_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}]}], 
   "}"}], "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"xA", ",", "xB", ",", "yA", ",", "yB"}], "}"}], "/.", 
   RowBox[{"FullSimplify", "[", 
    RowBox[{
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"u", "[", 
           RowBox[{
            RowBox[{"2", "\[Sigma]x"}], ",", "\[Alpha]", ",", 
            RowBox[{"Y", "+", "LyA", "+", "LyB"}], ",", 
            RowBox[{"pA", "+", "pB"}], ",", "z", ",", "1"}], "]"}], "-", 
          RowBox[{"u", "[", 
           RowBox[{"\[Sigma]x", ",", "\[Alpha]", ",", 
            RowBox[{"Y", "-", "yB", "+", "LyB"}], ",", "pB", ",", "z", ",", 
            RowBox[{"1", "-", "xA"}]}], "]"}]}], "==", "0"}], "&&", 
        RowBox[{
         RowBox[{
          RowBox[{"u", "[", 
           RowBox[{
            RowBox[{"2", "\[Sigma]x"}], ",", "\[Alpha]", ",", 
            RowBox[{"Y", "+", "LyA", "+", "LyB"}], ",", 
            RowBox[{"pA", "+", "pB"}], ",", "z", ",", "1"}], "]"}], "-", 
          RowBox[{"u", "[", 
           RowBox[{"\[Sigma]x", ",", "\[Alpha]", ",", 
            RowBox[{"yA", "+", "LyA"}], ",", "pA", ",", "z", ",", "xB"}], 
           "]"}]}], "==", "0"}], "&&", 
        RowBox[{
         RowBox[{
          RowBox[{"u", "[", 
           RowBox[{
            RowBox[{"2", "\[Sigma]y"}], ",", "\[Beta]", ",", 
            RowBox[{"1", "+", "LxA", "+", "LxB"}], ",", 
            RowBox[{"rA", "+", "rB"}], ",", "q", ",", "Y"}], "]"}], "-", 
          RowBox[{"u", "[", 
           RowBox[{"\[Sigma]y", ",", "\[Beta]", ",", 
            RowBox[{"1", "-", "xB", "+", "LxB"}], ",", "rB", ",", "q", ",", 
            RowBox[{"Y", "-", "yA"}]}], "]"}]}], "==", "0"}], "&&", 
        RowBox[{
         RowBox[{
          RowBox[{"u", "[", 
           RowBox[{
            RowBox[{"2", "\[Sigma]y"}], ",", "\[Beta]", ",", 
            RowBox[{"1", "+", "LxA", "+", "LxB"}], ",", 
            RowBox[{"rA", "+", "rB"}], ",", "q", ",", "Y"}], "]"}], "-", 
          RowBox[{"u", "[", 
           RowBox[{"\[Sigma]y", ",", "\[Beta]", ",", 
            RowBox[{"xA", "+", "LxA"}], ",", "rA", ",", "q", ",", "yB"}], 
           "]"}]}], "==", "0"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"xA", ",", "xB", ",", "yA", ",", "yB"}], "}"}]}], "]"}], "[", 
     RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.9620460642424297`*^9, 3.962046445172262*^9}, {
  3.962046483775606*^9, 3.962046484785091*^9}, {3.9620470752302628`*^9, 
  3.962047091362479*^9}, {3.962048783881861*^9, 3.962048824992587*^9}, {
  3.962048875191809*^9, 3.962048884671428*^9}, {3.962048935232052*^9, 
  3.962049007299182*^9}, {3.962049039939381*^9, 3.962049091471079*^9}, {
  3.962052138130747*^9, 3.962052200886145*^9}, {3.962145697631231*^9, 
  3.962145703861395*^9}},
 CellLabel->"In[12]:=",ExpressionUUID->"1da6c989-f96d-ae4b-a88d-e89c92f3d6b2"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   FractionBox[
    RowBox[{
     RowBox[{
      RowBox[{"-", "pA"}], " ", "q"}], "+", 
     RowBox[{"LyA", " ", "q", " ", "\[Alpha]"}], "+", 
     RowBox[{"rB", " ", "\[Alpha]"}], "+", 
     RowBox[{"q", " ", "Y", " ", "\[Alpha]"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]"}], "-", 
     RowBox[{"LxB", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
     RowBox[{"q", " ", "\[Sigma]x"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Sigma]y"}]}], 
    RowBox[{
     RowBox[{"q", " ", "z"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]], ",", 
   FractionBox[
    RowBox[{
     RowBox[{"pB", " ", "q"}], "+", 
     RowBox[{"q", " ", 
      RowBox[{"(", 
       RowBox[{"z", "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"LyB", "+", "Y"}], ")"}], " ", "\[Alpha]"}], "-", 
        "\[Sigma]x"}], ")"}]}], "+", 
     RowBox[{"\[Alpha]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "rA"}], "+", 
        RowBox[{"LxA", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
    RowBox[{
     RowBox[{"q", " ", "z"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]], ",", 
   FractionBox[
    RowBox[{
     RowBox[{
      RowBox[{"-", "rA"}], " ", "z"}], "+", 
     RowBox[{"pB", " ", "\[Beta]"}], "-", 
     RowBox[{"\[Beta]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"LyB", "+", "Y"}], ")"}], " ", "\[Alpha]"}], "+", 
        "\[Sigma]x"}], ")"}]}], "+", 
     RowBox[{"z", " ", 
      RowBox[{"(", 
       RowBox[{"\[Beta]", "+", 
        RowBox[{"LxA", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
    RowBox[{
     RowBox[{"q", " ", "z"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]], ",", 
   FractionBox[
    RowBox[{
     RowBox[{"rB", " ", "z"}], "+", 
     RowBox[{"q", " ", "Y", " ", "z"}], "-", 
     RowBox[{"\[Beta]", " ", 
      RowBox[{"(", 
       RowBox[{"pA", "+", "z", "+", 
        RowBox[{"LxB", " ", "z"}], "-", 
        RowBox[{"LyA", " ", "\[Alpha]"}], "-", "\[Sigma]x"}], ")"}]}], "-", 
     RowBox[{"z", " ", "\[Sigma]y"}]}], 
    RowBox[{
     RowBox[{"q", " ", "z"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]]}], "}"}]], "Output",
 CellChangeTimes->{3.9620463713827534`*^9, 3.962046449125364*^9, 
  3.962046488806045*^9, 3.962047096372818*^9, 3.962049097620228*^9, 
  3.9620522144610987`*^9, 3.962140084576605*^9, 3.962143012443809*^9, 
  3.96214532215801*^9, 3.962145569037853*^9, 3.962145713025776*^9, 
  3.962185837865505*^9, 3.962187010695334*^9, 3.962228686208305*^9, 
  3.962259235716057*^9, 3.9622602186886024`*^9, 3.962262923090419*^9, 
  3.962264075544269*^9, 3.962264568756148*^9, 3.962266030471817*^9, 
  3.9623039139868317`*^9, 3.962304274429085*^9, 3.962304845167874*^9, 
  3.962305247584996*^9, 3.9623055279402657`*^9, 3.9623056812574635`*^9, 
  3.962306054196059*^9, 3.9623068687108097`*^9, 3.962307093508865*^9, 
  3.962307186986271*^9, 3.96230778843767*^9, 3.9623079736917*^9, 
  3.962311148455462*^9, 3.962597419678249*^9, 3.962599332649212*^9, 
  3.962599952131092*^9, 3.9713410644957466`*^9},
 CellLabel->"Out[12]=",ExpressionUUID->"5d8a33c3-7c88-a247-bdf8-2ef28d1f7750"]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"xAhatD", "[", 
   RowBox[{
   "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_", 
    ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", "rB_",
     ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], ":=", 
  RowBox[{"xAhatSol", "[", 
   RowBox[{
   "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA", ",", 
    "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", "rB", ",", 
    "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"xBhatD", "[", 
   RowBox[{
   "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_", 
    ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", "rB_",
     ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], ":=", 
  RowBox[{"xBhatSol", "[", 
   RowBox[{
   "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA", ",", 
    "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", "rB", ",", 
    "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"yAhatD", "[", 
   RowBox[{
   "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_", 
    ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", "rB_",
     ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], ":=", 
  RowBox[{"yAhatSol", "[", 
   RowBox[{
   "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA", ",", 
    "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", "rB", ",", 
    "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"yBhatD", "[", 
   RowBox[{
   "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_", 
    ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", "rB_",
     ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], ":=", 
  RowBox[{"yBhatSol", "[", 
   RowBox[{
   "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA", ",", 
    "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", "rB", ",", 
    "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.9620470992914104`*^9, 3.962047135786175*^9}, {
  3.962047253044231*^9, 3.962047345823465*^9}, {3.96204911473761*^9, 
  3.962049213284173*^9}, {3.9620522254270267`*^9, 3.962052256529449*^9}, {
  3.9621457222533817`*^9, 3.962145759098678*^9}},
 CellLabel->"In[13]:=",ExpressionUUID->"8272f2e7-4d98-de46-9543-3c85c50de20a"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Profits", "Section",
 CellChangeTimes->{{3.9625926739086456`*^9, 
  3.962592676039238*^9}},ExpressionUUID->"58113e13-913c-294c-af9e-\
4f344f889b47"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"ProfitA", "[", 
   RowBox[{
   "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_", 
    ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", "rB_",
     ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], "=", 
  RowBox[{"FullSimplify", "[", 
   RowBox[{
    RowBox[{"pA", "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"xAhatD", "[", 
        RowBox[{
        "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA", 
         ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", "rB", 
         ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], "]"}], "+", 
       "LxA"}], ")"}]}], "+", 
    RowBox[{"rA", "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"yAhatD", "[", 
        RowBox[{
        "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA", 
         ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", "rB", 
         ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], "]"}], "+", 
       "LyA"}], ")"}]}]}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ProfitB", "[", 
   RowBox[{
   "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_", 
    ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", "rB_",
     ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], "=", 
  RowBox[{"FullSimplify", "[", 
   RowBox[{
    RowBox[{"pB", "*", 
     RowBox[{"(", 
      RowBox[{"1", "-", 
       RowBox[{"xBhatD", "[", 
        RowBox[{
        "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA", 
         ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", "rB", 
         ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], "]"}], "+", 
       "LxB"}], ")"}]}], "+", 
    RowBox[{"rB", "*", 
     RowBox[{"(", 
      RowBox[{"Y", "-", 
       RowBox[{"yBhatD", "[", 
        RowBox[{
        "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA", 
         ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", "rB", 
         ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], "]"}], "+", 
       "LyB"}], ")"}]}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.961798121398548*^9, 3.961798235780928*^9}, {
  3.961798342471131*^9, 3.961798360532295*^9}, {3.962045319098345*^9, 
  3.962045425044321*^9}, {3.9620466461385937`*^9, 3.962046653099682*^9}, {
  3.96204692880357*^9, 3.962047008259701*^9}, {3.962047359847229*^9, 
  3.9620474010080147`*^9}, {3.96204923132929*^9, 3.962049387348341*^9}, {
  3.962051282960911*^9, 3.962051295046182*^9}, {3.962052280683801*^9, 
  3.962052296728684*^9}, {3.962053007184748*^9, 3.962053028866721*^9}, {
  3.9620531127514534`*^9, 3.962053125328158*^9}, {3.962599367108364*^9, 
  3.962599378879145*^9}},
 CellLabel->"In[17]:=",ExpressionUUID->"287828de-3953-2148-ae4e-5eb78dbcdafa"],

Cell[BoxData[
 RowBox[{
  RowBox[{"LyA", " ", "rA"}], "+", 
  FractionBox[
   RowBox[{
    RowBox[{
     RowBox[{"-", 
      SuperscriptBox["pA", "2"]}], " ", "q"}], "+", 
    RowBox[{"pA", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"LxA", " ", "q", " ", "z"}], "+", 
       RowBox[{"LyA", " ", "q", " ", "\[Alpha]"}], "+", 
       RowBox[{"rB", " ", "\[Alpha]"}], "+", 
       RowBox[{"q", " ", "Y", " ", "\[Alpha]"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}], "-", 
       RowBox[{"LxA", " ", "\[Alpha]", " ", "\[Beta]"}], "-", 
       RowBox[{"LxB", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
       RowBox[{"q", " ", "\[Sigma]x"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Sigma]y"}]}], ")"}]}], "+", 
    RowBox[{"rA", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "rA"}], " ", "z"}], "+", 
       RowBox[{"pB", " ", "\[Beta]"}], "-", 
       RowBox[{"\[Beta]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"LyB", "+", "Y"}], ")"}], " ", "\[Alpha]"}], "+", 
          "\[Sigma]x"}], ")"}]}], "+", 
       RowBox[{"z", " ", 
        RowBox[{"(", 
         RowBox[{"\[Beta]", "+", 
          RowBox[{"LxA", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
      ")"}]}]}], 
   RowBox[{
    RowBox[{"q", " ", "z"}], "-", 
    RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]]}]], "Output",
 CellChangeTimes->{3.961798449054104*^9, 3.961969352350975*^9, 
  3.962043471814867*^9, 3.962045428086643*^9, 3.962047403321275*^9, 
  3.962049294633621*^9, 3.962049392496542*^9, 3.962051299585455*^9, 
  3.962052310316113*^9, 3.962053142786314*^9, 3.9621353220584755`*^9, 
  3.962137512010639*^9, 3.962592759502922*^9, 3.9625982962268105`*^9, 
  3.962598717140911*^9, 3.962598994350029*^9, 3.962599383533565*^9, 
  3.962599952892927*^9, 3.9713410654759407`*^9},
 CellLabel->"Out[17]=",ExpressionUUID->"eec37dfd-6fc2-c942-9f3b-2d3b474828a2"],

Cell[BoxData[
 RowBox[{
  RowBox[{"pB", " ", 
   RowBox[{"(", 
    RowBox[{"1", "+", "LxB", "-", 
     FractionBox[
      RowBox[{
       RowBox[{"pB", " ", "q"}], "+", 
       RowBox[{"q", " ", 
        RowBox[{"(", 
         RowBox[{"z", "-", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"LyB", "+", "Y"}], ")"}], " ", "\[Alpha]"}], "-", 
          "\[Sigma]x"}], ")"}]}], "+", 
       RowBox[{"\[Alpha]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "rA"}], "+", 
          RowBox[{"LxA", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
      RowBox[{
       RowBox[{"q", " ", "z"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]]}], ")"}]}], "+", 
  RowBox[{"rB", " ", 
   RowBox[{"(", 
    RowBox[{"LyB", "+", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", "rB"}], " ", "z"}], "+", 
       RowBox[{"pA", " ", "\[Beta]"}], "-", 
       RowBox[{"\[Beta]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"LyA", "+", "Y"}], ")"}], " ", "\[Alpha]"}], "+", 
          "\[Sigma]x"}], ")"}]}], "+", 
       RowBox[{"z", " ", 
        RowBox[{"(", 
         RowBox[{"\[Beta]", "+", 
          RowBox[{"LxB", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
      RowBox[{
       RowBox[{"q", " ", "z"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]]}], ")"}]}]}]], "Output",
 CellChangeTimes->{3.961798449054104*^9, 3.961969352350975*^9, 
  3.962043471814867*^9, 3.962045428086643*^9, 3.962047403321275*^9, 
  3.962049294633621*^9, 3.962049392496542*^9, 3.962051299585455*^9, 
  3.962052310316113*^9, 3.962053142786314*^9, 3.9621353220584755`*^9, 
  3.962137512010639*^9, 3.962592759502922*^9, 3.9625982962268105`*^9, 
  3.962598717140911*^9, 3.962598994350029*^9, 3.962599383533565*^9, 
  3.962599952892927*^9, 3.9713410661516075`*^9},
 CellLabel->"Out[18]=",ExpressionUUID->"122abe37-9b85-8949-a38d-2c7d7a00de03"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Profit Max->Best Responses", "Section",
 CellChangeTimes->{{3.9625926895258045`*^9, 
  3.962592701553026*^9}},ExpressionUUID->"e55b14d4-82f5-b349-828b-\
ce20e0dafd22"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"pABR", "[", 
      RowBox[{
      "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", 
       "pA_", ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", 
       ",", "rB_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], 
      "]"}], ",", 
     RowBox[{"rABR", "[", 
      RowBox[{
      "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", 
       "pA_", ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", 
       ",", "rB_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], 
      "]"}]}], "}"}], "=", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"pA", ",", "rA"}], "}"}], "/.", 
    RowBox[{"FullSimplify", "[", 
     RowBox[{
      RowBox[{"Solve", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"D", "[", 
           RowBox[{
            RowBox[{"ProfitA", "[", 
             RowBox[{
             "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
              "pA", ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", 
              ",", "rB", ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], 
             "]"}], ",", "pA"}], "]"}], "==", "0"}], "&&", 
         RowBox[{
          RowBox[{"D", "[", 
           RowBox[{
            RowBox[{"ProfitA", "[", 
             RowBox[{
             "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
              "pA", ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", 
              ",", "rB", ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], 
             "]"}], ",", "rA"}], "]"}], "==", "0"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"pA", ",", "rA"}], "}"}]}], "]"}], "[", 
      RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.9617985154102535`*^9, 3.961798536326399*^9}, {
  3.9617986028839664`*^9, 3.961798732366387*^9}, {3.961798766689533*^9, 
  3.961798771872376*^9}, {3.961798850802464*^9, 3.9617988966045036`*^9}, {
  3.961799278030489*^9, 3.9617993042877464`*^9}, {3.961799348851307*^9, 
  3.961799352180528*^9}, {3.961799382507313*^9, 3.961799389367586*^9}, {
  3.9617994271370316`*^9, 3.961799599998396*^9}, {3.9617996372333775`*^9, 
  3.961799664332033*^9}, {3.961969360706627*^9, 3.961969372737391*^9}, {
  3.962045444280233*^9, 3.962045499280096*^9}, {3.962047413827917*^9, 
  3.9620474327454243`*^9}, {3.962049451437627*^9, 3.962049523463964*^9}},
 CellLabel->"In[19]:=",ExpressionUUID->"5e8bcf18-5645-9d45-a88a-b5798abcbb00"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   FractionBox[
    RowBox[{
     RowBox[{"LxA", " ", "q", " ", "z"}], "+", 
     RowBox[{"LyA", " ", "q", " ", "\[Alpha]"}], "+", 
     RowBox[{"rB", " ", "\[Alpha]"}], "+", 
     RowBox[{"q", " ", "Y", " ", "\[Alpha]"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]"}], "-", 
     RowBox[{"LxA", " ", "\[Alpha]", " ", "\[Beta]"}], "-", 
     RowBox[{"LxB", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
     RowBox[{"q", " ", "\[Sigma]x"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Sigma]y"}]}], 
    RowBox[{"2", " ", "q"}]], ",", 
   FractionBox[
    RowBox[{
     RowBox[{"pB", " ", "\[Beta]"}], "+", 
     RowBox[{"LyA", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "z"}], "-", 
        RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "-", 
     RowBox[{"\[Beta]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"LyB", "+", "Y"}], ")"}], " ", "\[Alpha]"}], "+", 
        "\[Sigma]x"}], ")"}]}], "+", 
     RowBox[{"z", " ", 
      RowBox[{"(", 
       RowBox[{"\[Beta]", "+", 
        RowBox[{"LxA", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
    RowBox[{"2", " ", "z"}]]}], "}"}]], "Output",
 CellChangeTimes->{3.961799668554655*^9, 3.961969377026984*^9, 
  3.962045503156263*^9, 3.962047453401764*^9, 3.9620495318172965`*^9, 
  3.962051308854582*^9, 3.9620523244570737`*^9, 3.962053155683133*^9, 
  3.962135323704719*^9, 3.962137513520811*^9, 3.9625927602053432`*^9, 
  3.9625982969953804`*^9, 3.96259871788468*^9, 3.962599000252018*^9, 
  3.962599395244522*^9, 3.962599953641182*^9, 3.9713410662761497`*^9},
 CellLabel->"Out[19]=",ExpressionUUID->"df0cdac6-ae92-604a-ba64-7da497281222"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{"pBBR", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_",
       ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", 
      "rB_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], ",", 
    RowBox[{"rBBR", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_",
       ",", "pB_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", 
      "rB_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}]}], 
   "}"}], "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"pB", ",", "rB"}], "}"}], "/.", 
   RowBox[{"FullSimplify", "[", 
    RowBox[{
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{"ProfitB", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "pA", ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", 
             ",", "rB", ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], 
            "]"}], ",", "pB"}], "]"}], "==", "0"}], "&&", 
        RowBox[{
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{"ProfitB", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "pA", ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", 
             ",", "rB", ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], 
            "]"}], ",", "rB"}], "]"}], "==", "0"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"pB", ",", "rB"}], "}"}]}], "]"}], "[", 
     RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.962047465548353*^9, 3.9620476147570553`*^9}, {
  3.962049565953263*^9, 3.962049595384735*^9}},
 CellLabel->"In[20]:=",ExpressionUUID->"ffd973ad-bde9-6242-ab8c-b055aca7147c"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   FractionBox[
    RowBox[{
     RowBox[{"LxB", " ", "q", " ", "z"}], "+", 
     RowBox[{"LyB", " ", "q", " ", "\[Alpha]"}], "+", 
     RowBox[{"rA", " ", "\[Alpha]"}], "+", 
     RowBox[{"q", " ", "Y", " ", "\[Alpha]"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]"}], "-", 
     RowBox[{"LxA", " ", "\[Alpha]", " ", "\[Beta]"}], "-", 
     RowBox[{"LxB", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
     RowBox[{"q", " ", "\[Sigma]x"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Sigma]y"}]}], 
    RowBox[{"2", " ", "q"}]], ",", 
   FractionBox[
    RowBox[{
     RowBox[{"pA", " ", "\[Beta]"}], "+", 
     RowBox[{"LyB", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "z"}], "-", 
        RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "-", 
     RowBox[{"\[Beta]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"LyA", "+", "Y"}], ")"}], " ", "\[Alpha]"}], "+", 
        "\[Sigma]x"}], ")"}]}], "+", 
     RowBox[{"z", " ", 
      RowBox[{"(", 
       RowBox[{"\[Beta]", "+", 
        RowBox[{"LxB", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
    RowBox[{"2", " ", "z"}]]}], "}"}]], "Output",
 CellChangeTimes->{{3.9620475958111477`*^9, 3.962047616309826*^9}, 
   3.962049602142953*^9, 3.9620513136662364`*^9, 3.96205234543394*^9, 
   3.96205317416782*^9, 3.962135323775251*^9, 3.9621375135913296`*^9, 
   3.9625927602502956`*^9, 3.962598297050356*^9, 3.962598717930451*^9, 
   3.962599000438357*^9, 3.9625993953914394`*^9, 3.9625999536834717`*^9, 
   3.9713410663297787`*^9},
 CellLabel->"Out[20]=",ExpressionUUID->"3c59924c-fd64-9a41-82e6-b939844d7bf0"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Best Responses - Eqm", "Section",
 CellChangeTimes->{{3.96259270912916*^9, 
  3.962592722338978*^9}},ExpressionUUID->"15f2b005-cb3f-2241-a947-\
9897c46b1b9a"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{"pAeq", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", 
      "LyB_"}], "]"}], ",", 
    RowBox[{"rAeq", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", 
      "LyB_"}], "]"}], ",", 
    RowBox[{"pBeq", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", 
      "LyB_"}], "]"}], ",", 
    RowBox[{"rBeq", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", 
      "LyB_"}], "]"}]}], "}"}], "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"pA", ",", "rA", ",", "pB", ",", "rB"}], "}"}], "/.", 
   RowBox[{
    RowBox[{"FullSimplify", "[", 
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"pA", "==", 
         RowBox[{"pABR", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", 
           "rB", ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], "]"}]}], "&&", 
        RowBox[{"rA", "==", 
         RowBox[{"rABR", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", 
           "rB", ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], "]"}]}], "&&", 
        RowBox[{"pB", "==", 
         RowBox[{"pBBR", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", 
           "rB", ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], "]"}]}], "&&", 
        RowBox[{"rB", "==", 
         RowBox[{"rBBR", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", 
           "rB", ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], "]"}]}]}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"pA", ",", "rA", ",", "pB", ",", "rB"}], "}"}]}], "]"}], 
     "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.9621353578535004`*^9, 3.962135378964207*^9}},
 CellLabel->"In[25]:=",ExpressionUUID->"5aacb4e4-5e09-b64d-b3f2-c7f1e9d05f33"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", 
    FractionBox[
     RowBox[{
      RowBox[{
       RowBox[{"-", "LyB"}], " ", "q", " ", "z", " ", "\[Alpha]"}], "-", 
      RowBox[{"2", " ", "q", " ", "Y", " ", "z", " ", "\[Alpha]"}], "+", 
      RowBox[{"z", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
      RowBox[{"LxB", " ", "z", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
      RowBox[{"LyB", " ", 
       SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]"}], "+", 
      RowBox[{"Y", " ", 
       SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]"}], "+", 
      RowBox[{"LyA", " ", "\[Alpha]", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "2"}], " ", "q", " ", "z"}], "+", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "+", 
      RowBox[{"2", " ", "LxA", " ", "z", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "q"}], " ", "z"}], "+", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "-", 
      RowBox[{"2", " ", "q", " ", "z", " ", "\[Sigma]x"}], "+", 
      RowBox[{"\[Alpha]", " ", "\[Beta]", " ", "\[Sigma]x"}], "+", 
      RowBox[{"z", " ", "\[Alpha]", " ", "\[Sigma]y"}]}], 
     RowBox[{
      RowBox[{"4", " ", "q", " ", "z"}], "-", 
      RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]]}], ",", 
   FractionBox[
    RowBox[{
     RowBox[{"2", " ", "LyA", " ", "q", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "z"}], "-", 
        RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1", "+", "LxA", "+", "LxB"}], ")"}], " ", "\[Beta]"}], "+",
         "\[Sigma]y"}], ")"}]}], "+", 
     RowBox[{"q", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"2", "+", 
           RowBox[{"2", " ", "LxA"}], "+", "LxB"}], ")"}], " ", "z", " ", 
         "\[Beta]"}], "-", 
        RowBox[{"\[Beta]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"LyB", "+", "Y"}], ")"}], " ", "\[Alpha]"}], "+", 
           "\[Sigma]x"}], ")"}]}], "+", 
        RowBox[{"2", " ", "z", " ", "\[Sigma]y"}]}], ")"}]}]}], 
    RowBox[{
     RowBox[{"4", " ", "q", " ", "z"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]], ",", 
   RowBox[{"-", 
    FractionBox[
     RowBox[{
      RowBox[{
       RowBox[{"-", "2"}], " ", "LyB", " ", "q", " ", "z", " ", "\[Alpha]"}], 
      "-", 
      RowBox[{"2", " ", "q", " ", "Y", " ", "z", " ", "\[Alpha]"}], "+", 
      RowBox[{"z", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
      RowBox[{"LxA", " ", "z", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
      RowBox[{"LyB", " ", 
       SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]"}], "+", 
      RowBox[{"Y", " ", 
       SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]"}], "+", 
      RowBox[{"2", " ", "LxB", " ", "z", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "q"}], " ", "z"}], "+", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "+", 
      RowBox[{"LyA", " ", "\[Alpha]", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "q"}], " ", "z"}], "+", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "-", 
      RowBox[{"2", " ", "q", " ", "z", " ", "\[Sigma]x"}], "+", 
      RowBox[{"\[Alpha]", " ", "\[Beta]", " ", "\[Sigma]x"}], "+", 
      RowBox[{"z", " ", "\[Alpha]", " ", "\[Sigma]y"}]}], 
     RowBox[{
      RowBox[{"4", " ", "q", " ", "z"}], "-", 
      RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]]}], ",", 
   FractionBox[
    RowBox[{
     RowBox[{"2", " ", "LyB", " ", "q", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "z"}], "-", 
        RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1", "+", "LxA", "+", "LxB"}], ")"}], " ", "\[Beta]"}], "+",
         "\[Sigma]y"}], ")"}]}], "+", 
     RowBox[{"q", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"2", "+", "LxA", "+", 
           RowBox[{"2", " ", "LxB"}]}], ")"}], " ", "z", " ", "\[Beta]"}], 
        "-", 
        RowBox[{"\[Beta]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"LyA", "+", "Y"}], ")"}], " ", "\[Alpha]"}], "+", 
           "\[Sigma]x"}], ")"}]}], "+", 
        RowBox[{"2", " ", "z", " ", "\[Sigma]y"}]}], ")"}]}]}], 
    RowBox[{
     RowBox[{"4", " ", "q", " ", "z"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]]}], "}"}]], "Output",
 CellChangeTimes->{3.9621353810989857`*^9, 3.962137515049118*^9, 
  3.962592760819527*^9, 3.96259829761878*^9, 3.9625990067293415`*^9, 
  3.96259940553084*^9, 3.9625999542582603`*^9, 3.9713410795175896`*^9},
 CellLabel->"Out[25]=",ExpressionUUID->"b3f21212-24bd-984a-9afd-849e10c687ac"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["If LxA=L, LxB=0, LyA=0, LyB=0", "Section",
 CellChangeTimes->{{3.971341232482464*^9, 
  3.971341271734215*^9}},ExpressionUUID->"eaf531c3-bdfe-1849-a57a-\
f69925b30ae8"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{"pAeqL", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "L_"}], "]"}], ",", 
    RowBox[{"rAeqL", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "L_"}], "]"}], ",", 
    RowBox[{"pBeqL", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "L_"}], "]"}], ",", 
    RowBox[{"rBeqL", "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "L_"}], "]"}]}], "}"}], "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"pA", ",", "rA", ",", "pB", ",", "rB"}], "}"}], "/.", 
   RowBox[{
    RowBox[{"FullSimplify", "[", 
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"pA", "==", 
         RowBox[{"pABR", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", 
           "rB", ",", "L", ",", "0", ",", "0", ",", "0"}], "]"}]}], "&&", 
        RowBox[{"rA", "==", 
         RowBox[{"rABR", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", 
           "rB", ",", "L", ",", "0", ",", "0", ",", "0"}], "]"}]}], "&&", 
        RowBox[{"pB", "==", 
         RowBox[{"pBBR", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", 
           "rB", ",", "L", ",", "0", ",", "0", ",", "0"}], "]"}]}], "&&", 
        RowBox[{"rB", "==", 
         RowBox[{"rBBR", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", 
           "rB", ",", "L", ",", "0", ",", "0", ",", "0"}], "]"}]}]}], ",", 
       RowBox[{"{", 
        RowBox[{"pA", ",", "rA", ",", "pB", ",", "rB"}], "}"}]}], "]"}], 
     "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.971341274837904*^9, 3.9713413130299606`*^9}},
 CellLabel->"In[26]:=",ExpressionUUID->"426004d3-426f-744b-8e72-dabc00fae52d"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   FractionBox[
    RowBox[{
     RowBox[{"2", " ", "L", " ", "z", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "z"}], "-", 
        RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "+", 
     RowBox[{"2", " ", "q", " ", "z", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "-", 
     RowBox[{"\[Alpha]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"\[Beta]", " ", 
         RowBox[{"(", 
          RowBox[{"z", "+", 
           RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "+", 
        RowBox[{"z", " ", "\[Sigma]y"}]}], ")"}]}]}], 
    RowBox[{
     RowBox[{"4", " ", "q", " ", "z"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]], ",", 
   RowBox[{"-", 
    FractionBox[
     RowBox[{
      RowBox[{"q", " ", "\[Beta]", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "-", 
      RowBox[{"2", " ", "q", " ", "z", " ", 
       RowBox[{"(", 
        RowBox[{"\[Beta]", "+", 
         RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}], "+", 
      RowBox[{"\[Alpha]", " ", "\[Beta]", " ", 
       RowBox[{"(", 
        RowBox[{"\[Beta]", "+", 
         RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
     RowBox[{
      RowBox[{"4", " ", "q", " ", "z"}], "-", 
      RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]]}], ",", 
   FractionBox[
    RowBox[{
     RowBox[{"2", " ", "q", " ", "z", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "-", 
     RowBox[{"\[Alpha]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"\[Beta]", " ", 
         RowBox[{"(", 
          RowBox[{"z", "+", 
           RowBox[{"L", " ", "z"}], "+", 
           RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "+", 
        RowBox[{"z", " ", "\[Sigma]y"}]}], ")"}]}]}], 
    RowBox[{
     RowBox[{"4", " ", "q", " ", "z"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]], ",", 
   FractionBox[
    RowBox[{
     RowBox[{
      RowBox[{"-", "\[Alpha]"}], " ", "\[Beta]", " ", 
      RowBox[{"(", 
       RowBox[{"\[Beta]", "+", 
        RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}], "+", 
     RowBox[{"q", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"2", "+", "L"}], ")"}], " ", "z", " ", "\[Beta]"}], "-", 
        RowBox[{"\[Beta]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "+", 
        RowBox[{"2", " ", "z", " ", "\[Sigma]y"}]}], ")"}]}]}], 
    RowBox[{
     RowBox[{"4", " ", "q", " ", "z"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]]}], "}"}]], "Output",
 CellChangeTimes->{3.9713413188781834`*^9},
 CellLabel->"Out[26]=",ExpressionUUID->"c35f5ddc-1bae-7240-a7b9-eb4ab2a0d4aa"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Specific Value in Basic Model", "Section",
 CellChangeTimes->{{3.9625927301459465`*^9, 3.96259274101165*^9}, {
  3.962598266671261*^9, 3.962598284964426*^9}, {3.9713411312178783`*^9, 
  3.9713411364592113`*^9}},ExpressionUUID->"19fb375e-efa7-834e-97ec-\
10cfe7b7a2ba"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"pAeq", "[", 
  RowBox[{
  "1.6", ",", "1.6", " ", ",", "0.6", ",", "1", ",", "1.4", ",", "1.6", ",", 
   "1.5", ",", "1", ",", "0", ",", "0", ",", "0"}], "]"}]], "Input",
 CellChangeTimes->{{3.96213539555825*^9, 3.9621353964157467`*^9}, {
   3.9621354400949306`*^9, 3.962135440958554*^9}, {3.962137543168427*^9, 
   3.962137548775348*^9}, 3.97134133558827*^9},
 CellLabel->"In[28]:=",ExpressionUUID->"5815f8ef-544d-264f-b39e-3bb52976e630"],

Cell[BoxData["1.1196029776674938`"], "Output",
 CellChangeTimes->{3.9621353984478664`*^9, 3.962135444148519*^9, 
  3.962137515084648*^9, 3.962137550892624*^9, 3.962592760839524*^9, 
  3.962598297652855*^9, 3.962599011231653*^9, 3.962599409107113*^9, 
  3.9625999542883224`*^9, 3.971341362739893*^9},
 CellLabel->"Out[28]=",ExpressionUUID->"12850f40-705a-fc42-952d-fd60d00b9927"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"rAeq", "[", 
  RowBox[{
  "1.6", ",", "1.6", " ", ",", "0.6", ",", "1", ",", "1.4", ",", "1.6", ",", 
   "1.5", ",", "1", ",", "0", ",", "0", ",", "0"}], "]"}]], "Input",
 CellChangeTimes->{{3.9621354498727093`*^9, 3.962135457402569*^9}, {
   3.962137554026596*^9, 3.962137557412948*^9}, 3.9713413383583736`*^9},
 CellLabel->"In[29]:=",ExpressionUUID->"7a63be0e-5932-7041-8d96-665aefd17e2a"],

Cell[BoxData["1.388089330024814`"], "Output",
 CellChangeTimes->{3.9621354587487946`*^9, 3.962137515092804*^9, 
  3.962137558613823*^9, 3.962592760865677*^9, 3.9625982976796913`*^9, 
  3.962599011290039*^9, 3.962599409167353*^9, 3.9625999543238125`*^9, 
  3.9713413666137905`*^9},
 CellLabel->"Out[29]=",ExpressionUUID->"04e313ed-031f-ed46-ad6e-5bc79eeb9947"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"pBeq", "[", 
  RowBox[{
  "1.6", ",", "1.6", " ", ",", "0.6", ",", "1", ",", "1.4", ",", "1.6", ",", 
   "1.5", ",", "1", ",", "0", ",", "0", ",", "0"}], "]"}]], "Input",
 CellChangeTimes->{{3.962135516666197*^9, 3.962135522456789*^9}, {
   3.962135553782789*^9, 3.9621356011168785`*^9}, {3.9621375003162746`*^9, 
   3.962137524634976*^9}, {3.962137561696306*^9, 3.962137565823191*^9}, 
   3.971341340630804*^9},
 CellLabel->"In[30]:=",ExpressionUUID->"37f57791-5eed-014d-8c2d-5b961029eeae"],

Cell[BoxData["0.49776674937965265`"], "Output",
 CellChangeTimes->{3.962135602129563*^9, 3.96213752764826*^9, 
  3.962137567729761*^9, 3.962592760891785*^9, 3.9625982976897144`*^9, 
  3.96259901129702*^9, 3.962599409197354*^9, 3.96259995435594*^9, 
  3.9713413681819515`*^9},
 CellLabel->"Out[30]=",ExpressionUUID->"a1666f0a-40da-cc41-908d-acdb756f1f9b"]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Find \[Sigma]x Upper Bound", "Chapter",
 CellChangeTimes->{{3.962599786716076*^9, 
  3.9625997943177624`*^9}},ExpressionUUID->"31228355-b033-5941-8aa4-\
15146b7ed2f5"],

Cell[CellGroupData[{

Cell["Best Responses - Only Eqm pB,rB", "Section",
 CellChangeTimes->{{3.962593025913101*^9, 3.9625930404085846`*^9}, {
  3.9625968095236664`*^9, 3.962596814231344*^9}, {3.96259848127075*^9, 
  3.9625984992756634`*^9}},ExpressionUUID->"be53c7ca-f7a6-394e-8ddd-\
26bff516b8ca"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{
     RowBox[{"pBeq", "'"}], "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", 
      "LyB_"}], "]"}], ",", 
    RowBox[{
     RowBox[{"rBeq", "'"}], "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "LxA_", ",", "LxB_", ",", "LyA_", ",", 
      "LyB_"}], "]"}]}], "}"}], "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"pB", ",", "rB"}], "}"}], "/.", 
   RowBox[{
    RowBox[{"FullSimplify", "[", 
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"pB", "==", 
         RowBox[{"pBBR", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", 
           "rB", ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], "]"}]}], "&&", 
        RowBox[{"rB", "==", 
         RowBox[{"rBBR", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", "pB", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", 
           "rB", ",", "LxA", ",", "LxB", ",", "LyA", ",", "LyB"}], "]"}]}]}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"pB", ",", "rB"}], "}"}]}], "]"}], "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.962229678406582*^9, 3.962229685041395*^9}, {
  3.962229745999991*^9, 3.962229752418083*^9}, {3.962229812211225*^9, 
  3.962229824086521*^9}, {3.962305006585934*^9, 3.962305009453964*^9}, {
  3.962307128248134*^9, 3.962307135784389*^9}, {3.962598512035223*^9, 
  3.962598513026663*^9}},
 CellLabel->"In[31]:=",ExpressionUUID->"f0cbbfc3-6237-8246-b486-bd41e4e20838"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   FractionBox[
    RowBox[{
     RowBox[{"LxB", " ", "q", " ", "z"}], "+", 
     RowBox[{"LyB", " ", "q", " ", "\[Alpha]"}], "+", 
     RowBox[{"rA", " ", "\[Alpha]"}], "+", 
     RowBox[{"q", " ", "Y", " ", "\[Alpha]"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Beta]"}], "-", 
     RowBox[{"LxA", " ", "\[Alpha]", " ", "\[Beta]"}], "-", 
     RowBox[{"LxB", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
     RowBox[{"q", " ", "\[Sigma]x"}], "-", 
     RowBox[{"\[Alpha]", " ", "\[Sigma]y"}]}], 
    RowBox[{"2", " ", "q"}]], ",", 
   FractionBox[
    RowBox[{
     RowBox[{"pA", " ", "\[Beta]"}], "+", 
     RowBox[{"LyB", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "z"}], "-", 
        RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "-", 
     RowBox[{"\[Beta]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"LyA", "+", "Y"}], ")"}], " ", "\[Alpha]"}], "+", 
        "\[Sigma]x"}], ")"}]}], "+", 
     RowBox[{"z", " ", 
      RowBox[{"(", 
       RowBox[{"\[Beta]", "+", 
        RowBox[{"LxB", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
    RowBox[{"2", " ", "z"}]]}], "}"}]], "Output",
 CellChangeTimes->{3.9622298270405903`*^9, 3.962259235989969*^9, 
  3.962260219053602*^9, 3.962262923457726*^9, 3.962264075884149*^9, 
  3.962264569057516*^9, 3.96226603080108*^9, 3.9623039143994904`*^9, 
  3.9623042750697517`*^9, 3.9623048455492706`*^9, 3.962305011395172*^9, 
  3.962305248134533*^9, 3.962305528371788*^9, 3.962305727701191*^9, 
  3.962306055789406*^9, 3.9623068701036873`*^9, 3.96230713974654*^9, 
  3.9623072296513405`*^9, 3.9623077899784355`*^9, 3.962307975194801*^9, 
  3.962311150517628*^9, 3.9625974205054398`*^9, 3.962599029523428*^9, 
  3.962599412428009*^9, 3.962599954390188*^9, 3.971341413387533*^9},
 CellLabel->"Out[31]=",ExpressionUUID->"53d859d5-bfcc-ad44-88cf-62f8d5753584"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["New Indifferent Consumers", "Section",
 CellChangeTimes->{{3.962317077406445*^9, 3.962317085484066*^9}, {
  3.9625930639261055`*^9, 3.9625930641797085`*^9}, {3.9625968201532*^9, 
  3.9625968403537693`*^9}},ExpressionUUID->"cd991f6d-b2a6-5646-b66c-\
1efca7415d69"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{
     RowBox[{"xAhat", "'"}], "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "pA_", ",", "Y_",
       ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", "LxA_", ",", 
      "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], ",", 
    RowBox[{
     RowBox[{"xBhat", "'"}], "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "pA_", ",", "Y_",
       ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", "LxA_", ",", 
      "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], ",", 
    RowBox[{
     RowBox[{"yAhat", "'"}], "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "pA_", ",", "Y_",
       ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", "LxA_", ",", 
      "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], ",", 
    RowBox[{
     RowBox[{"yBhat", "'"}], "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "pA_", ",", "Y_",
       ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", "LxA_", ",", 
      "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}]}], "}"}], "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"xA", ",", "xB", ",", "yA", ",", "yB"}], "}"}], "/.", 
   RowBox[{
    RowBox[{"FullSimplify", "[", 
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"xA", "==", 
         RowBox[{"xAhatD", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", 
           RowBox[{
            RowBox[{"pBeq", "'"}], "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "LxA", ",", "LxB", ",", 
             "LyA", ",", "LyB"}], "]"}], ",", "z", ",", "q", ",", "\[Beta]", 
           ",", "rA", ",", 
           RowBox[{
            RowBox[{"rBeq", "'"}], "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "LxA", ",", "LxB", ",", 
             "LyA", ",", "LyB"}], "]"}], ",", "L", ",", "0", ",", "0", ",", 
           "0"}], "]"}]}], "&&", 
        RowBox[{"xB", "==", 
         RowBox[{"xBhatD", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", 
           RowBox[{
            RowBox[{"pBeq", "'"}], "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "LxA", ",", "LxB", ",", 
             "LyA", ",", "LyB"}], "]"}], ",", "z", ",", "q", ",", "\[Beta]", 
           ",", "rA", ",", 
           RowBox[{
            RowBox[{"rBeq", "'"}], "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "LxA", ",", "LxB", ",", 
             "LyA", ",", "LyB"}], "]"}], ",", "L", ",", "0", ",", "0", ",", 
           "0"}], "]"}]}], "&&", 
        RowBox[{"yA", "==", 
         RowBox[{"yAhatD", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", 
           RowBox[{
            RowBox[{"pBeq", "'"}], "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "LxA", ",", "LxB", ",", 
             "LyA", ",", "LyB"}], "]"}], ",", "z", ",", "q", ",", "\[Beta]", 
           ",", "rA", ",", 
           RowBox[{
            RowBox[{"rBeq", "'"}], "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "LxA", ",", "LxB", ",", 
             "LyA", ",", "LyB"}], "]"}], ",", "L", ",", "0", ",", "0", ",", 
           "0"}], "]"}]}], "&&", 
        RowBox[{"yB", "==", 
         RowBox[{"yBhatD", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA",
            ",", 
           RowBox[{
            RowBox[{"pBeq", "'"}], "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "LxA", ",", "LxB", ",", 
             "LyA", ",", "LyB"}], "]"}], ",", "z", ",", "q", ",", "\[Beta]", 
           ",", "rA", ",", 
           RowBox[{
            RowBox[{"rBeq", "'"}], "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "LxA", ",", "LxB", ",", 
             "LyA", ",", "LyB"}], "]"}], ",", "L", ",", "0", ",", "0", ",", 
           "0"}], "]"}]}]}], ",", 
       RowBox[{"{", 
        RowBox[{"xA", ",", "xB", ",", "yA", ",", "yB"}], "}"}]}], "]"}], 
     "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.962229108797825*^9, 3.962229235696377*^9}, {
  3.962229603478674*^9, 3.9622296328322277`*^9}, {3.9622298443939915`*^9, 
  3.962229905016941*^9}, {3.962229936184551*^9, 3.962229959628748*^9}, {
  3.9622300226597843`*^9, 3.962230029905474*^9}, {3.962305035445135*^9, 
  3.962305057350384*^9}, {3.962307242319599*^9, 3.962307251795439*^9}, {
  3.9623072898418674`*^9, 3.962307291817171*^9}},
 CellLabel->"In[32]:=",ExpressionUUID->"64219419-dd83-4e49-862c-54ab2a83ce12"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", 
    FractionBox[
     RowBox[{
      RowBox[{"2", " ", "pA", " ", "q", " ", "z"}], "-", 
      RowBox[{"LyB", " ", "q", " ", "z", " ", "\[Alpha]"}], "-", 
      RowBox[{"2", " ", "q", " ", "Y", " ", "z", " ", "\[Alpha]"}], "-", 
      RowBox[{"pA", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
      RowBox[{"z", " ", "\[Alpha]", " ", "\[Beta]"}], "-", 
      RowBox[{"LxB", " ", "z", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
      RowBox[{"LyA", " ", 
       SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]"}], "+", 
      RowBox[{"LyB", " ", 
       SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]"}], "+", 
      RowBox[{"Y", " ", 
       SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]"}], "-", 
      RowBox[{"2", " ", "q", " ", "z", " ", "\[Sigma]x"}], "+", 
      RowBox[{"\[Alpha]", " ", "\[Beta]", " ", "\[Sigma]x"}], "+", 
      RowBox[{"z", " ", "\[Alpha]", " ", "\[Sigma]y"}]}], 
     RowBox[{
      RowBox[{"2", " ", "q", " ", 
       SuperscriptBox["z", "2"]}], "-", 
      RowBox[{"2", " ", "z", " ", "\[Alpha]", " ", "\[Beta]"}]}]]}], ",", 
   FractionBox[
    RowBox[{
     RowBox[{"q", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"2", "+", "LxB"}], ")"}], " ", "z"}], "+", 
        RowBox[{"LyB", " ", "\[Alpha]"}], "-", 
        RowBox[{"Y", " ", "\[Alpha]"}], "-", "\[Sigma]x"}], ")"}]}], "-", 
     RowBox[{"\[Alpha]", " ", 
      RowBox[{"(", 
       RowBox[{"rA", "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1", "-", 
           RowBox[{"2", " ", "L"}], "+", "LxA", "+", "LxB"}], ")"}], " ", 
         "\[Beta]"}], "-", "\[Sigma]y"}], ")"}]}]}], 
    RowBox[{
     RowBox[{"2", " ", "q", " ", "z"}], "-", 
     RowBox[{"2", " ", "\[Alpha]", " ", "\[Beta]"}]}]], ",", 
   FractionBox[
    RowBox[{
     RowBox[{
      RowBox[{"-", "\[Alpha]"}], " ", "\[Beta]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "rA"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1", "+", "LxA", "+", "LxB"}], ")"}], " ", "\[Beta]"}], "+",
         "\[Sigma]y"}], ")"}]}], "+", 
     RowBox[{"q", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"-", "2"}], " ", "rA", " ", "z"}], "+", 
        RowBox[{"\[Beta]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"2", "+", 
              RowBox[{"2", " ", "L"}], "+", "LxB"}], ")"}], " ", "z"}], "+", 
           RowBox[{"LyB", " ", "\[Alpha]"}], "-", 
           RowBox[{"Y", " ", "\[Alpha]"}], "-", "\[Sigma]x"}], ")"}]}], "+", 
        RowBox[{"2", " ", "z", " ", "\[Sigma]y"}]}], ")"}]}]}], 
    RowBox[{"2", " ", "q", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"q", " ", "z"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]], ",", 
   FractionBox[
    RowBox[{
     RowBox[{"2", " ", "q", " ", "Y", " ", "z"}], "+", 
     RowBox[{"LyB", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "z"}], "-", 
        RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "-", 
     RowBox[{"\[Beta]", " ", 
      RowBox[{"(", 
       RowBox[{"pA", "+", "z", "-", 
        RowBox[{"LxB", " ", "z"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"LyA", "+", "Y"}], ")"}], " ", "\[Alpha]"}], "-", 
        "\[Sigma]x"}], ")"}]}], "-", 
     RowBox[{"z", " ", "\[Sigma]y"}]}], 
    RowBox[{
     RowBox[{"2", " ", "q", " ", "z"}], "-", 
     RowBox[{"2", " ", "\[Alpha]", " ", "\[Beta]"}]}]]}], "}"}]], "Output",
 CellChangeTimes->{3.9622299118834305`*^9, 3.962229964389139*^9, 
  3.962230033174879*^9, 3.962259236328117*^9, 3.962260219384124*^9, 
  3.9622629238470364`*^9, 3.9622640761755753`*^9, 3.962264569355945*^9, 
  3.962266031139766*^9, 3.962303914801832*^9, 3.9623042756651096`*^9, 
  3.962304846014979*^9, 3.9623050622499046`*^9, 3.962305248637512*^9, 
  3.962305528907376*^9, 3.962305741139969*^9, 3.962306056566599*^9, 
  3.9623068710470695`*^9, 3.962307301009985*^9, 3.9623077910641994`*^9, 
  3.962307976167015*^9, 3.962311151446129*^9, 3.962597420956932*^9, 
  3.9625994284001713`*^9, 3.962599954788378*^9, 3.9713414137985516`*^9},
 CellLabel->"Out[32]=",ExpressionUUID->"a551a93a-f316-c548-b74d-cf357c984c9e"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["New ProfitA", "Section",
 CellChangeTimes->{{3.9623170974196644`*^9, 3.962317121549764*^9}, {
  3.9623171665624046`*^9, 3.962317169838688*^9}, {3.962596848315098*^9, 
  3.962596852618889*^9}},ExpressionUUID->"195018d1-e59b-624f-996a-\
cebc149953f4"],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"ProfitA", "=", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"xAhat", "'"}], "+", "LxA"}], ")"}], "*", "pA"}], "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"yAhat", "'"}], "+", "LyA"}], ")"}], "*", "rA"}]}]}], 
   ")"}]}]], "Input",
 CellChangeTimes->{{3.9622301291024437`*^9, 3.962230212362852*^9}, {
  3.962230251956529*^9, 3.962230252891188*^9}},
 EmphasizeSyntaxErrors->True,
 CellLabel->"In[33]:=",ExpressionUUID->"d13025f6-611f-8a40-a1bd-90d159f1cb41"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"ProfitA", "'"}], "[", 
   RowBox[{
   "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_", 
    ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "rA_", ",", "LxA_", ",", 
    "LxB_", ",", "LyA_", ",", "LyB_"}], "]"}], "=", 
  RowBox[{"FullSimplify", "[", 
   RowBox[{
    RowBox[{"pA", "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"xAhat", "'"}], "[", 
        RowBox[{
        "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "pA", ",", "Y", 
         ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", "LxA", ",", 
         "LxB", ",", "LyA", ",", "LyB"}], "]"}], "+", "LxA"}], ")"}]}], "+", 
    RowBox[{"rA", "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"yAhat", "'"}], "[", 
        RowBox[{
        "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "pA", ",", "Y", 
         ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", "LxA", ",", 
         "LxB", ",", "LyA", ",", "LyB"}], "]"}], "+", "LyA"}], ")"}]}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.962230292674731*^9, 3.962230431790966*^9}, 
   3.9622593062003784`*^9, {3.962305754241974*^9, 3.962305780949932*^9}},
 CellLabel->"In[33]:=",ExpressionUUID->"ca6399b8-d9f7-cc43-9c6c-5580c076bcc3"],

Cell[BoxData[
 RowBox[{
  RowBox[{"pA", " ", 
   RowBox[{"(", 
    RowBox[{"LxA", "-", 
     FractionBox[
      RowBox[{
       RowBox[{"2", " ", "pA", " ", "q", " ", "z"}], "-", 
       RowBox[{"LyB", " ", "q", " ", "z", " ", "\[Alpha]"}], "-", 
       RowBox[{"2", " ", "q", " ", "Y", " ", "z", " ", "\[Alpha]"}], "-", 
       RowBox[{"pA", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
       RowBox[{"z", " ", "\[Alpha]", " ", "\[Beta]"}], "-", 
       RowBox[{"LxB", " ", "z", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
       RowBox[{"LyA", " ", 
        SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]"}], "+", 
       RowBox[{"LyB", " ", 
        SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]"}], "+", 
       RowBox[{"Y", " ", 
        SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]"}], "-", 
       RowBox[{"2", " ", "q", " ", "z", " ", "\[Sigma]x"}], "+", 
       RowBox[{"\[Alpha]", " ", "\[Beta]", " ", "\[Sigma]x"}], "+", 
       RowBox[{"z", " ", "\[Alpha]", " ", "\[Sigma]y"}]}], 
      RowBox[{
       RowBox[{"2", " ", "q", " ", 
        SuperscriptBox["z", "2"]}], "-", 
       RowBox[{"2", " ", "z", " ", "\[Alpha]", " ", "\[Beta]"}]}]]}], ")"}]}],
   "+", 
  RowBox[{"rA", " ", 
   RowBox[{"(", 
    RowBox[{"LyA", "+", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", "\[Alpha]"}], " ", "\[Beta]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "rA"}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"1", "+", "LxA", "+", "LxB"}], ")"}], " ", "\[Beta]"}], 
          "+", "\[Sigma]y"}], ")"}]}], "+", 
       RowBox[{"q", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"-", "2"}], " ", "rA", " ", "z"}], "+", 
          RowBox[{"\[Beta]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{"2", "+", 
                RowBox[{"2", " ", "L"}], "+", "LxB"}], ")"}], " ", "z"}], "+", 
             RowBox[{"LyB", " ", "\[Alpha]"}], "-", 
             RowBox[{"Y", " ", "\[Alpha]"}], "-", "\[Sigma]x"}], ")"}]}], "+", 
          RowBox[{"2", " ", "z", " ", "\[Sigma]y"}]}], ")"}]}]}], 
      RowBox[{"2", " ", "q", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"q", " ", "z"}], "-", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]]}], 
    ")"}]}]}]], "Output",
 CellChangeTimes->{3.962230444053444*^9, 3.9622592367107525`*^9, 
  3.962259315298258*^9, 3.962260219686235*^9, 3.96226292425896*^9, 
  3.9622640764556713`*^9, 3.9622645696396503`*^9, 3.96226603150024*^9, 
  3.9623039153178864`*^9, 3.962304276185255*^9, 3.962304846448363*^9, 
  3.9623050757279854`*^9, 3.962305249379553*^9, 3.9623055295878563`*^9, 
  3.9623057857321873`*^9, 3.9623060591292534`*^9, 3.962306873310236*^9, 
  3.962307328329497*^9, 3.9623077930905247`*^9, 3.962307978123537*^9, 
  3.9623111540561905`*^9, 3.9625974221889095`*^9, 3.962599491398195*^9, 
  3.962599956042364*^9, 3.9713414149794884`*^9},
 CellLabel->"Out[33]=",ExpressionUUID->"c1b0a2f9-a635-6d47-9409-07fa31e3a155"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["New rA", "Section",
 CellChangeTimes->{{3.9625968658181114`*^9, 
  3.9625968686900616`*^9}},ExpressionUUID->"e8996b59-8377-1745-b19d-\
5cbad812efca"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{"rABR", "'"}], "[", 
    RowBox[{
    "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "pA_", 
     ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "LxA_", ",", "LxB_", ",", 
     "LyA_", ",", "LyB_"}], "]"}], "}"}], "=", 
  RowBox[{
   RowBox[{"{", "rA", "}"}], "/.", 
   RowBox[{"FullSimplify", "[", 
    RowBox[{
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"ProfitA", "'"}], "[", 
           RowBox[{
           "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
            "pA", ",", "z", ",", "q", ",", "\[Beta]", ",", "rA", ",", "LxA", 
            ",", "LxB", ",", "LyA", ",", "LyB"}], "]"}], ",", "rA"}], "]"}], "==",
         "0"}], ",", 
       RowBox[{"{", "rA", "}"}]}], "]"}], "[", 
     RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.962230475397278*^9, 3.962230527756212*^9}, 
   3.962230889759186*^9, 3.962259311731571*^9, 3.962259615689407*^9, 
   3.9623050834596977`*^9},
 CellLabel->"In[34]:=",ExpressionUUID->"56c43505-5841-714b-a949-3ba28d803d6a"],

Cell[BoxData[
 RowBox[{"{", 
  FractionBox[
   RowBox[{
    RowBox[{"2", " ", "LyA", " ", "q", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"q", " ", "z"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "+", 
    RowBox[{"q", " ", "\[Beta]", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"2", "+", 
          RowBox[{"2", " ", "L"}], "+", "LxB"}], ")"}], " ", "z"}], "+", 
       RowBox[{"LyB", " ", "\[Alpha]"}], "-", 
       RowBox[{"Y", " ", "\[Alpha]"}], "-", "\[Sigma]x"}], ")"}]}], "+", 
    RowBox[{"2", " ", "q", " ", "z", " ", "\[Sigma]y"}], "-", 
    RowBox[{"\[Alpha]", " ", "\[Beta]", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "+", "LxA", "+", "LxB"}], ")"}], " ", "\[Beta]"}], "+", 
       "\[Sigma]y"}], ")"}]}]}], 
   RowBox[{
    RowBox[{"4", " ", "q", " ", "z"}], "-", 
    RowBox[{"2", " ", "\[Alpha]", " ", "\[Beta]"}]}]], "}"}]], "Output",
 CellChangeTimes->{3.962230592378126*^9, 3.962230893760931*^9, 
  3.962259236850424*^9, 3.962259321306886*^9, 3.962259618535433*^9, 
  3.9622602198272057`*^9, 3.962262924463383*^9, 3.9622640765730553`*^9, 
  3.96226456975194*^9, 3.9622660326518784`*^9, 3.962303915651771*^9, 
  3.962304276459183*^9, 3.962304846633227*^9, 3.962305086507789*^9, 
  3.962305249598425*^9, 3.962305529823297*^9, 3.962305798875918*^9, 
  3.962306059350115*^9, 3.9623068735336304`*^9, 3.96230733533971*^9, 
  3.9623077932829285`*^9, 3.962307978331555*^9, 3.9623111543718014`*^9, 
  3.9625974223436184`*^9, 3.962599501457926*^9, 3.962599956214985*^9, 
  3.97134141512142*^9},
 CellLabel->"Out[34]=",ExpressionUUID->"3de31723-4264-7748-88bf-f9b189369f47"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["New ProfitA<=Original ProfitA -> New pA Upper Bound ", "Section",
 CellChangeTimes->{{3.9623172355808735`*^9, 3.962317256642473*^9}, {
  3.962597093752102*^9, 3.9625971299026966`*^9}, {3.962597164514887*^9, 
  3.962597196782667*^9}, {3.962597292700323*^9, 
  3.962597325581412*^9}},ExpressionUUID->"09599678-18b0-0446-ae13-\
db4e325c73bc"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"pA", "'"}], "[", 
   RowBox[{
   "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
    ",", "q_", ",", "\[Beta]_", ",", "L_"}], "]"}], "=", 
  RowBox[{"pA", "/.", 
   RowBox[{"FullSimplify", "[", 
    RowBox[{
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"ProfitA", "'"}], "[", 
         RowBox[{
         "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "pA", 
          ",", "z", ",", "q", ",", "\[Beta]", ",", 
          RowBox[{
           RowBox[{"rABR", "'"}], "[", 
           RowBox[{
           "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
            "pA", ",", "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",",
             "0", ",", "0"}], "]"}], ",", "L", ",", "0", ",", "0", ",", "0"}],
          "]"}], "==", 
        RowBox[{"ProfitA", "[", 
         RowBox[{
         "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
          RowBox[{"pAeq", "[", 
           RowBox[{
           "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z",
             ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",", 
            "0"}], "]"}], ",", 
          RowBox[{"pBeq", "[", 
           RowBox[{
           "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z",
             ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",", 
            "0"}], "]"}], ",", "z", ",", "q", ",", "\[Beta]", ",", 
          RowBox[{"rAeq", "[", 
           RowBox[{
           "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z",
             ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",", 
            "0"}], "]"}], ",", 
          RowBox[{"rBeq", "[", 
           RowBox[{
           "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z",
             ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",", 
            "0"}], "]"}], ",", "L", ",", "0", ",", "0", ",", "0"}], "]"}]}], 
       ",", "pA"}], "]"}], "[", 
     RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.962307694527355*^9, 3.962307701899277*^9}},
 CellLabel->"In[35]:=",ExpressionUUID->"31f1adab-b60d-ac44-9856-6da97ad627db"],

Cell[BoxData[
 RowBox[{"-", 
  RowBox[{
   FractionBox["1", 
    RowBox[{
     RowBox[{"4", " ", "q", " ", "z"}], "-", 
     RowBox[{"2", " ", "\[Alpha]", " ", "\[Beta]"}]}]], 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "2"}], " ", "q", " ", "z", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"L", " ", "z"}], "+", 
        RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "+", 
     RowBox[{"\[Alpha]", " ", "\[Beta]", " ", 
      RowBox[{"(", 
       RowBox[{"z", "+", 
        RowBox[{"2", " ", "L", " ", "z"}], "+", 
        RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "+", 
     RowBox[{"z", " ", "\[Alpha]", " ", "\[Sigma]y"}], "-", 
     RowBox[{"q", " ", 
      SuperscriptBox["z", "2"], " ", 
      RowBox[{"\[Sqrt]", 
       RowBox[{"(", 
        RowBox[{
         FractionBox["1", 
          RowBox[{"q", " ", 
           SuperscriptBox["z", "2"], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{"q", " ", "z"}], "-", 
              RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], "2"], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "4"}], " ", "q", " ", "z"}], "+", 
              RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], "2"]}]], 
         RowBox[{
          SuperscriptBox["\[Alpha]", "2"], " ", 
          SuperscriptBox["\[Beta]", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             SuperscriptBox["L", "2"], " ", "z", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"4", " ", 
                SuperscriptBox["q", "3"], " ", 
                SuperscriptBox["z", "3"]}], "+", 
               RowBox[{"4", " ", "q", " ", "z", " ", "\[Alpha]", " ", 
                RowBox[{"(", 
                 RowBox[{"\[Alpha]", "-", "\[Beta]"}], ")"}], " ", 
                SuperscriptBox["\[Beta]", "2"]}], "+", 
               RowBox[{
                SuperscriptBox["\[Alpha]", "2"], " ", 
                SuperscriptBox["\[Beta]", "4"]}], "+", 
               RowBox[{"4", " ", 
                SuperscriptBox["q", "2"], " ", 
                SuperscriptBox["z", "2"], " ", "\[Beta]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"-", "2"}], " ", "\[Alpha]"}], "+", "\[Beta]"}], 
                 ")"}]}]}], ")"}]}], "+", 
            RowBox[{"4", " ", 
             SuperscriptBox["q", "3"], " ", 
             SuperscriptBox["z", "2"], " ", 
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{
                RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}], 
              "2"]}], "+", 
            RowBox[{"z", " ", 
             SuperscriptBox["\[Alpha]", "2"], " ", 
             SuperscriptBox["\[Beta]", "2"], " ", 
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}], "+", 
            RowBox[{"q", " ", "\[Alpha]", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"\[Alpha]", " ", 
                SuperscriptBox["\[Beta]", "2"], " ", 
                SuperscriptBox[
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}], 
                 "2"]}], "+", 
               RowBox[{"2", " ", "z", " ", "\[Beta]", " ", 
                RowBox[{"(", 
                 RowBox[{"\[Alpha]", "+", "\[Beta]"}], ")"}], " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}], 
                " ", 
                RowBox[{"(", 
                 RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "+", 
               RowBox[{
                SuperscriptBox["z", "2"], " ", 
                RowBox[{"(", 
                 RowBox[{"\[Alpha]", "-", 
                  RowBox[{"4", " ", "\[Beta]"}]}], ")"}], " ", 
                SuperscriptBox[
                 RowBox[{"(", 
                  RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}]}], 
              ")"}]}], "+", 
            RowBox[{"2", " ", "L", " ", "z", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"4", " ", 
                SuperscriptBox["q", "3"], " ", 
                SuperscriptBox["z", "2"], " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}],
                "+", 
               RowBox[{
                SuperscriptBox["\[Alpha]", "2"], " ", 
                SuperscriptBox["\[Beta]", "3"], " ", 
                RowBox[{"(", 
                 RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "-", 
               RowBox[{"2", " ", 
                SuperscriptBox["q", "2"], " ", "z", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"\[Beta]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"3", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}], 
                   " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], 
                    ")"}]}], "+", 
                  RowBox[{"z", " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Alpha]", "-", 
                    RowBox[{"2", " ", "\[Beta]"}]}], ")"}], " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], ")"}]}],
                "+", 
               RowBox[{"q", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"\[Beta]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"2", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}], 
                   " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], 
                    ")"}]}], "+", 
                  RowBox[{"2", " ", "z", " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Alpha]", "-", 
                    RowBox[{"2", " ", "\[Beta]"}]}], ")"}], " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], 
                 ")"}]}]}], ")"}]}], "-", 
            RowBox[{
             SuperscriptBox["q", "2"], " ", "z", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{
                SuperscriptBox["Y", "2"], " ", 
                SuperscriptBox["\[Alpha]", "2"], " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"4", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                " ", "\[Beta]"}], "+", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"4", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                " ", "\[Beta]", " ", 
                SuperscriptBox["\[Sigma]x", "2"]}], "+", 
               RowBox[{"4", " ", "z", " ", 
                RowBox[{"(", 
                 RowBox[{"\[Alpha]", "+", "\[Beta]"}], ")"}], " ", 
                "\[Sigma]x", " ", 
                RowBox[{"(", 
                 RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "-", 
               RowBox[{"4", " ", 
                SuperscriptBox["z", "2"], " ", 
                SuperscriptBox[
                 RowBox[{"(", 
                  RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}], "+", 
               RowBox[{"2", " ", "Y", " ", "\[Alpha]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"4", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                   " ", "\[Beta]", " ", "\[Sigma]x"}], "+", 
                  RowBox[{"2", " ", "z", " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Alpha]", "+", "\[Beta]"}], ")"}], " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], 
                 ")"}]}]}], ")"}]}]}], ")"}]}]}], ")"}]}]}], "+", 
     RowBox[{"z", " ", "\[Alpha]", " ", "\[Beta]", " ", 
      RowBox[{"\[Sqrt]", 
       RowBox[{"(", 
        RowBox[{
         FractionBox["1", 
          RowBox[{"q", " ", 
           SuperscriptBox["z", "2"], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{"q", " ", "z"}], "-", 
              RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], "2"], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "4"}], " ", "q", " ", "z"}], "+", 
              RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], "2"]}]], 
         RowBox[{
          SuperscriptBox["\[Alpha]", "2"], " ", 
          SuperscriptBox["\[Beta]", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             SuperscriptBox["L", "2"], " ", "z", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"4", " ", 
                SuperscriptBox["q", "3"], " ", 
                SuperscriptBox["z", "3"]}], "+", 
               RowBox[{"4", " ", "q", " ", "z", " ", "\[Alpha]", " ", 
                RowBox[{"(", 
                 RowBox[{"\[Alpha]", "-", "\[Beta]"}], ")"}], " ", 
                SuperscriptBox["\[Beta]", "2"]}], "+", 
               RowBox[{
                SuperscriptBox["\[Alpha]", "2"], " ", 
                SuperscriptBox["\[Beta]", "4"]}], "+", 
               RowBox[{"4", " ", 
                SuperscriptBox["q", "2"], " ", 
                SuperscriptBox["z", "2"], " ", "\[Beta]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"-", "2"}], " ", "\[Alpha]"}], "+", "\[Beta]"}], 
                 ")"}]}]}], ")"}]}], "+", 
            RowBox[{"4", " ", 
             SuperscriptBox["q", "3"], " ", 
             SuperscriptBox["z", "2"], " ", 
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{
                RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}], 
              "2"]}], "+", 
            RowBox[{"z", " ", 
             SuperscriptBox["\[Alpha]", "2"], " ", 
             SuperscriptBox["\[Beta]", "2"], " ", 
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}], "+", 
            RowBox[{"q", " ", "\[Alpha]", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"\[Alpha]", " ", 
                SuperscriptBox["\[Beta]", "2"], " ", 
                SuperscriptBox[
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}], 
                 "2"]}], "+", 
               RowBox[{"2", " ", "z", " ", "\[Beta]", " ", 
                RowBox[{"(", 
                 RowBox[{"\[Alpha]", "+", "\[Beta]"}], ")"}], " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}], 
                " ", 
                RowBox[{"(", 
                 RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "+", 
               RowBox[{
                SuperscriptBox["z", "2"], " ", 
                RowBox[{"(", 
                 RowBox[{"\[Alpha]", "-", 
                  RowBox[{"4", " ", "\[Beta]"}]}], ")"}], " ", 
                SuperscriptBox[
                 RowBox[{"(", 
                  RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}]}], 
              ")"}]}], "+", 
            RowBox[{"2", " ", "L", " ", "z", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"4", " ", 
                SuperscriptBox["q", "3"], " ", 
                SuperscriptBox["z", "2"], " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}],
                "+", 
               RowBox[{
                SuperscriptBox["\[Alpha]", "2"], " ", 
                SuperscriptBox["\[Beta]", "3"], " ", 
                RowBox[{"(", 
                 RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "-", 
               RowBox[{"2", " ", 
                SuperscriptBox["q", "2"], " ", "z", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"\[Beta]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"3", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}], 
                   " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], 
                    ")"}]}], "+", 
                  RowBox[{"z", " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Alpha]", "-", 
                    RowBox[{"2", " ", "\[Beta]"}]}], ")"}], " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], ")"}]}],
                "+", 
               RowBox[{"q", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"\[Beta]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"2", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}], 
                   " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], 
                    ")"}]}], "+", 
                  RowBox[{"2", " ", "z", " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Alpha]", "-", 
                    RowBox[{"2", " ", "\[Beta]"}]}], ")"}], " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], 
                 ")"}]}]}], ")"}]}], "-", 
            RowBox[{
             SuperscriptBox["q", "2"], " ", "z", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{
                SuperscriptBox["Y", "2"], " ", 
                SuperscriptBox["\[Alpha]", "2"], " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"4", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                " ", "\[Beta]"}], "+", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"4", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                " ", "\[Beta]", " ", 
                SuperscriptBox["\[Sigma]x", "2"]}], "+", 
               RowBox[{"4", " ", "z", " ", 
                RowBox[{"(", 
                 RowBox[{"\[Alpha]", "+", "\[Beta]"}], ")"}], " ", 
                "\[Sigma]x", " ", 
                RowBox[{"(", 
                 RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "-", 
               RowBox[{"4", " ", 
                SuperscriptBox["z", "2"], " ", 
                SuperscriptBox[
                 RowBox[{"(", 
                  RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}], "+", 
               RowBox[{"2", " ", "Y", " ", "\[Alpha]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"4", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                   " ", "\[Beta]", " ", "\[Sigma]x"}], "+", 
                  RowBox[{"2", " ", "z", " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Alpha]", "+", "\[Beta]"}], ")"}], " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], 
                 ")"}]}]}], ")"}]}]}], ")"}]}]}], ")"}]}]}]}], 
    ")"}]}]}]], "Output",
 CellChangeTimes->{3.962307735158624*^9, 3.962307824114441*^9, 
  3.962308001925621*^9, 3.962311180451098*^9, 3.9625974356012745`*^9, 
  3.962599537846998*^9, 3.9625999677382507`*^9, 3.9713414268964233`*^9},
 CellLabel->"Out[35]=",ExpressionUUID->"90cc936a-c9b5-6344-b64e-68c625c9e60d"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\[Sigma]x=pA -> \[Sigma]x Upper Bound", "Section",
 CellChangeTimes->{{3.9623173250328484`*^9, 3.962317328522843*^9}, {
  3.962597354558596*^9, 
  3.9625973866232204`*^9}},ExpressionUUID->"affa3178-3bb7-7b44-b932-\
16da0c4dbf2b"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Sigma]x", "'"}], "[", 
   RowBox[{
   "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", ",", "q_", ",", 
    "\[Beta]_", ",", "L_"}], "]"}], "=", 
  RowBox[{"\[Sigma]x", "/.", 
   RowBox[{"FullSimplify", "[", 
    RowBox[{
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{"\[Sigma]x", "==", 
        RowBox[{
         RowBox[{"pA", "'"}], "[", 
         RowBox[{
         "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z", 
          ",", "q", ",", "\[Beta]", ",", "L"}], "]"}]}], ",", "\[Sigma]x"}], 
      "]"}], "[", 
     RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.962314833172943*^9, 3.9623149151321793`*^9}},
 CellLabel->"In[37]:=",ExpressionUUID->"ec5783c5-fa28-334f-83f0-c5e0f7ccfcbf"],

Cell[BoxData[
 TemplateBox[{
  "Solve", "nongen", 
   "\"\:53ef\:80fd\:5b58\:5728\:4f7f\:67d0\:4e9b\:89e3\:6216\:6240\:6709\:89e3\
\:90fd\:65e0\:6548\:7684\:53c2\:6570\:503c.\"", 2, 37, 3, 
   17181794852937535281, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.9713414644627666`*^9},
 CellLabel->
  "\:6b63\:5728\:8ba1\:7b97In[37]:=",ExpressionUUID->"567307af-5a62-7b49-9e73-\
9b84b3b850f5"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"64", " ", 
     SuperscriptBox["q", "4"], " ", 
     SuperscriptBox["z", "4"], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"L", " ", "z"}], "+", 
       RowBox[{"Y", " ", "\[Alpha]"}]}], ")"}]}], "+", 
    RowBox[{"8", " ", 
     SuperscriptBox["q", "2"], " ", 
     SuperscriptBox["z", "2"], " ", 
     SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"4", "+", 
          RowBox[{"11", " ", "L"}]}], ")"}], " ", "z", " ", "\[Beta]"}], "+", 
       RowBox[{"7", " ", "Y", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
       RowBox[{"4", " ", "z", " ", "\[Sigma]y"}]}], ")"}]}], "+", 
    RowBox[{
     SuperscriptBox["\[Alpha]", "3"], " ", 
     SuperscriptBox["\[Beta]", "3"], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Beta]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"2", " ", "\[Alpha]", " ", 
           RowBox[{"(", 
            RowBox[{"z", "+", 
             RowBox[{"2", " ", "L", " ", "z"}], "+", 
             RowBox[{"Y", " ", "\[Alpha]"}]}], ")"}]}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"1", "+", "L"}], ")"}], " ", "z", " ", "\[Beta]"}]}], 
         ")"}]}], "+", 
       RowBox[{"z", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"2", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}], " ", 
        "\[Sigma]y"}]}], ")"}]}], "-", 
    RowBox[{"32", " ", 
     SuperscriptBox["q", "3"], " ", 
     SuperscriptBox["z", "3"], " ", "\[Alpha]", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"3", " ", "Y", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
       RowBox[{"z", " ", 
        RowBox[{"(", 
         RowBox[{"\[Beta]", "+", 
          RowBox[{"4", " ", "L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], 
         ")"}]}]}], ")"}]}], "-", 
    RowBox[{"q", " ", "z", " ", 
     SuperscriptBox["\[Alpha]", "2"], " ", 
     SuperscriptBox["\[Beta]", "2"], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"Y", " ", "\[Alpha]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"16", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], " ", 
        "\[Beta]"}], "+", 
       RowBox[{"2", " ", "z", " ", "\[Beta]", " ", 
        RowBox[{"(", 
         RowBox[{"\[Beta]", "+", 
          RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}], "+", 
       RowBox[{"4", " ", "z", " ", "\[Alpha]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"3", "+", 
             RowBox[{"7", " ", "L"}]}], ")"}], " ", "\[Beta]"}], "+", 
          RowBox[{"3", " ", "\[Sigma]y"}]}], ")"}]}]}], ")"}]}], "-", 
    RowBox[{"2", " ", "q", " ", "z", " ", 
     RowBox[{"\[Sqrt]", 
      RowBox[{"(", 
       RowBox[{
        FractionBox["1", 
         RowBox[{
          SuperscriptBox["q", "2"], " ", 
          SuperscriptBox["z", "2"]}]], 
        RowBox[{
         SuperscriptBox["\[Alpha]", "2"], " ", 
         SuperscriptBox["\[Beta]", "2"], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", " ", "q", " ", "z"}], "-", 
           RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", " ", 
            SuperscriptBox["L", "2"], " ", 
            SuperscriptBox["z", "2"], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"64", " ", 
               SuperscriptBox["q", "5"], " ", 
               SuperscriptBox["z", "5"]}], "+", 
              RowBox[{"4", " ", 
               SuperscriptBox["q", "3"], " ", 
               SuperscriptBox["z", "3"], " ", "\[Alpha]", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"53", " ", "\[Alpha]"}], "-", 
                 RowBox[{"2", " ", "\[Beta]"}]}], ")"}], " ", 
               SuperscriptBox["\[Beta]", "2"]}], "+", 
              RowBox[{"q", " ", "z", " ", 
               SuperscriptBox["\[Alpha]", "3"], " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"24", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
               " ", 
               SuperscriptBox["\[Beta]", "4"]}], "-", 
              RowBox[{"2", " ", 
               SuperscriptBox["\[Alpha]", "5"], " ", 
               SuperscriptBox["\[Beta]", "5"]}], "+", 
              RowBox[{"4", " ", 
               SuperscriptBox["q", "4"], " ", 
               SuperscriptBox["z", "4"], " ", "\[Beta]", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "48"}], " ", "\[Alpha]"}], "+", "\[Beta]"}], 
                ")"}]}], "+", 
              RowBox[{"2", " ", 
               SuperscriptBox["q", "2"], " ", 
               SuperscriptBox["z", "2"], " ", 
               SuperscriptBox["\[Alpha]", "2"], " ", 
               SuperscriptBox["\[Beta]", "3"], " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "53"}], " ", "\[Alpha]"}], "+", 
                 RowBox[{"3", " ", "\[Beta]"}]}], ")"}]}]}], ")"}]}], "+", 
           RowBox[{"L", " ", "z", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"4", " ", "q", " ", "z"}], "-", 
              RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"64", " ", 
               SuperscriptBox["q", "4"], " ", "Y", " ", 
               SuperscriptBox["z", "4"], " ", "\[Alpha]"}], "+", 
              RowBox[{"2", " ", 
               SuperscriptBox["q", "2"], " ", 
               SuperscriptBox["z", "2"], " ", "\[Alpha]", " ", "\[Beta]", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"Y", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"56", " ", "\[Alpha]"}], "+", 
                    RowBox[{"5", " ", "\[Beta]"}]}], ")"}]}], "+", 
                 RowBox[{"2", " ", "z", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"14", " ", "\[Alpha]"}], "-", 
                    RowBox[{"3", " ", "\[Beta]"}]}], ")"}], " ", 
                  RowBox[{"(", 
                   RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], ")"}]}], 
              "-", 
              RowBox[{"2", " ", "q", " ", "z", " ", 
               SuperscriptBox["\[Alpha]", "2"], " ", 
               SuperscriptBox["\[Beta]", "2"], " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"3", " ", "Y", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"6", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}]}],
                  "+", 
                 RowBox[{"2", " ", "z", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"7", " ", "\[Alpha]"}], "-", 
                    RowBox[{"2", " ", "\[Beta]"}]}], ")"}], " ", 
                  RowBox[{"(", 
                   RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], ")"}]}], 
              "+", 
              RowBox[{
               SuperscriptBox["\[Alpha]", "3"], " ", 
               SuperscriptBox["\[Beta]", "3"], " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"Y", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"4", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}]}],
                  "+", 
                 RowBox[{"z", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"4", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                  " ", 
                  RowBox[{"(", 
                   RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], ")"}]}], 
              "-", 
              RowBox[{"8", " ", 
               SuperscriptBox["q", "3"], " ", 
               SuperscriptBox["z", "3"], " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"Y", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"18", " ", "\[Alpha]"}], "+", "\[Beta]"}], 
                   ")"}]}], "+", 
                 RowBox[{"z", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"4", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                  " ", 
                  RowBox[{"(", 
                   RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], 
                ")"}]}]}], ")"}]}], "+", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{
                RowBox[{"-", "4"}], " ", "q", " ", "z"}], "+", 
               RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], "2"], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"8", " ", 
               SuperscriptBox["q", "3"], " ", 
               SuperscriptBox["Y", "2"], " ", 
               SuperscriptBox["z", "3"], " ", 
               SuperscriptBox["\[Alpha]", "2"]}], "+", 
              RowBox[{"2", " ", 
               SuperscriptBox["q", "2"], " ", 
               SuperscriptBox["z", "2"], " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  SuperscriptBox["Y", "2"], " ", 
                  SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{
                    RowBox[{"-", "6"}], " ", "\[Alpha]"}], "+", "\[Beta]"}], 
                   ")"}]}], "-", 
                 RowBox[{"2", " ", "Y", " ", "z", " ", "\[Alpha]", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"2", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}], 
                  " ", 
                  RowBox[{"(", 
                   RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "+", 
                 RowBox[{
                  SuperscriptBox["z", "2"], " ", 
                  SuperscriptBox[
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}]}], 
                ")"}]}], "-", 
              RowBox[{
               SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  SuperscriptBox["Y", "2"], " ", 
                  SuperscriptBox["\[Alpha]", "3"], " ", 
                  SuperscriptBox["\[Beta]", "2"]}], "+", 
                 RowBox[{"Y", " ", "z", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"2", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}], 
                  " ", 
                  RowBox[{"(", 
                   RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "+", 
                 RowBox[{
                  SuperscriptBox["z", "2"], " ", 
                  RowBox[{"(", 
                   RowBox[{"\[Alpha]", "-", "\[Beta]"}], ")"}], " ", 
                  SuperscriptBox[
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}]}], 
                ")"}]}], "+", 
              RowBox[{"q", " ", "z", " ", "\[Alpha]", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  SuperscriptBox["Y", "2"], " ", 
                  SuperscriptBox["\[Alpha]", "2"], " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"6", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                  " ", 
                  SuperscriptBox["\[Beta]", "2"]}], "+", 
                 RowBox[{"Y", " ", "z", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"8", " ", "\[Alpha]"}], "+", 
                    RowBox[{"3", " ", "\[Beta]"}]}], ")"}], " ", 
                  RowBox[{"(", 
                   RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "+", 
                 RowBox[{"2", " ", 
                  SuperscriptBox["z", "2"], " ", 
                  RowBox[{"(", 
                   RowBox[{"\[Alpha]", "-", "\[Beta]"}], ")"}], " ", 
                  SuperscriptBox[
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}]}], 
                ")"}]}]}], ")"}]}]}], ")"}]}]}], ")"}]}]}]}], ")"}], "/", 
  RowBox[{"(", 
   RowBox[{"q", " ", "z", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"64", " ", 
       SuperscriptBox["q", "3"], " ", 
       SuperscriptBox["z", "3"]}], "-", 
      RowBox[{"96", " ", 
       SuperscriptBox["q", "2"], " ", 
       SuperscriptBox["z", "2"], " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
      RowBox[{"48", " ", "q", " ", "z", " ", 
       SuperscriptBox["\[Alpha]", "2"], " ", 
       SuperscriptBox["\[Beta]", "2"]}], "-", 
      RowBox[{
       SuperscriptBox["\[Alpha]", "2"], " ", 
       SuperscriptBox["\[Beta]", "3"], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"8", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}]}]}], ")"}]}],
    ")"}]}]], "Output",
 CellChangeTimes->{3.962314948060044*^9, 3.9625974494196644`*^9, 
  3.9625995647056484`*^9, 3.962599980798897*^9, 3.9713414749233418`*^9},
 CellLabel->"Out[37]=",ExpressionUUID->"6bdae327-aa66-2c49-83cf-a9b64c6e5cfa"]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Find \[Sigma]x Lower Bound", "Chapter",
 CellChangeTimes->{{3.962599822476206*^9, 
  3.9625998301268406`*^9}},ExpressionUUID->"1a2d8b0b-902e-e44e-ab23-\
ff4488f72778"],

Cell[CellGroupData[{

Cell["\<\
Indifferent Consumers from Original Eqm, only include \[Sigma]x,\[Sigma]y,\
\[Alpha],\[Beta],z,q,Y,L\
\>", "Section",
 CellChangeTimes->{{3.962324618431689*^9, 3.9623246408977637`*^9}, {
  3.962597660374436*^9, 3.962597731874563*^9}, {3.962597802350254*^9, 
  3.962597803136095*^9}, {3.9625998460262585`*^9, 
  3.9625998770925617`*^9}},ExpressionUUID->"304c1d60-8da9-c14b-b4f8-\
da5fe23ed6ba"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{
     RowBox[{"xAhat", "''"}], "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "L_"}], "]"}], ",", 
    RowBox[{
     RowBox[{"xBhat", "''"}], "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "L_"}], "]"}], ",", 
    RowBox[{
     RowBox[{"yAhat", "''"}], "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "L_"}], "]"}], ",", 
    RowBox[{
     RowBox[{"yBhat", "''"}], "[", 
     RowBox[{
     "\[Sigma]x_", ",", "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", 
      ",", "q_", ",", "\[Beta]_", ",", "L_"}], "]"}]}], "}"}], "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"xA", ",", "xB", ",", "yA", ",", "yB"}], "}"}], "/.", 
   RowBox[{"FullSimplify", "[", 
    RowBox[{
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"xA", "==", 
         RowBox[{"xAhatD", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
           RowBox[{"pAeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", 
           RowBox[{"pBeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", "z", ",", "q", ",", "\[Beta]", ",", 
           RowBox[{"rAeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", 
           RowBox[{"rBeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", "L", ",", "0", ",", "0", ",", "0"}], "]"}]}],
         "&&", 
        RowBox[{"xB", "==", 
         RowBox[{"xBhatD", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
           RowBox[{"pAeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", 
           RowBox[{"pBeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", "z", ",", "q", ",", "\[Beta]", ",", 
           RowBox[{"rAeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", 
           RowBox[{"rBeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", "L", ",", "0", ",", "0", ",", "0"}], "]"}]}],
         "&&", 
        RowBox[{"yA", "==", 
         RowBox[{"yAhatD", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
           RowBox[{"pAeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", 
           RowBox[{"pBeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", "z", ",", "q", ",", "\[Beta]", ",", 
           RowBox[{"rAeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", 
           RowBox[{"rBeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", "L", ",", "0", ",", "0", ",", "0"}], "]"}]}],
         "&&", 
        RowBox[{"yB", "==", 
         RowBox[{"yBhatD", "[", 
          RowBox[{
          "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
           RowBox[{"pAeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", 
           RowBox[{"pBeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", "z", ",", "q", ",", "\[Beta]", ",", 
           RowBox[{"rAeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", 
           RowBox[{"rBeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", "L", ",", "0", ",", "0", ",", "0"}], 
          "]"}]}]}], ",", 
       RowBox[{"{", 
        RowBox[{"xA", ",", "xB", ",", "yA", ",", "yB"}], "}"}]}], "]"}], "[", 
     RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.962322735334509*^9, 3.962322941331934*^9}, {
  3.962322982880131*^9, 3.96232310848588*^9}, {3.962323652055401*^9, 
  3.9623237117000446`*^9}, {3.9625997007752748`*^9, 3.962599706697171*^9}},
 CellLabel->"In[38]:=",ExpressionUUID->"3ace2da8-68f0-8d42-a9ae-1f9137c2cb87"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   FractionBox[
    RowBox[{
     RowBox[{
      RowBox[{"-", "L"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"q", " ", "z"}], "-", 
        RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"2", " ", "q", " ", "z"}], "-", 
        RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "+", 
     RowBox[{"q", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"2", " ", "q", " ", "z", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "-", 
        RowBox[{"\[Alpha]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"\[Beta]", " ", 
            RowBox[{"(", 
             RowBox[{"z", "+", 
              RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], 
           "+", 
           RowBox[{"z", " ", "\[Sigma]y"}]}], ")"}]}]}], ")"}]}]}], 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"q", " ", "z"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"4", " ", "q", " ", "z"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]], ",", 
   FractionBox[
    RowBox[{
     RowBox[{
      SuperscriptBox["\[Alpha]", "2"], " ", 
      SuperscriptBox["\[Beta]", "2"]}], "+", 
     RowBox[{"2", " ", 
      SuperscriptBox["q", "2"], " ", "z", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"2", " ", "z"}], "-", 
        RowBox[{"Y", " ", "\[Alpha]"}], "-", "\[Sigma]x"}], ")"}]}], "+", 
     RowBox[{"q", " ", "\[Alpha]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"\[Beta]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "4"}], "+", "L"}], ")"}], " ", "z"}], "+", 
           RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "+", 
        RowBox[{"z", " ", "\[Sigma]y"}]}], ")"}]}]}], 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"q", " ", "z"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"4", " ", "q", " ", "z"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]], ",", 
   FractionBox[
    RowBox[{"z", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "q"}], " ", "\[Beta]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "+", 
       RowBox[{"2", " ", "q", " ", "z", " ", 
        RowBox[{"(", 
         RowBox[{"\[Beta]", "+", 
          RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]", " ", 
        RowBox[{"(", 
         RowBox[{"\[Beta]", "+", 
          RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
      ")"}]}], 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"q", " ", "z"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"4", " ", "q", " ", "z"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]], ",", 
   FractionBox[
    RowBox[{
     RowBox[{"4", " ", 
      SuperscriptBox["q", "2"], " ", "Y", " ", 
      SuperscriptBox["z", "2"]}], "-", 
     RowBox[{"q", " ", "z", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"\[Beta]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"2", "+", "L"}], ")"}], " ", "z"}], "+", 
           RowBox[{"4", " ", "Y", " ", "\[Alpha]"}], "-", "\[Sigma]x"}], 
          ")"}]}], "+", 
        RowBox[{"2", " ", "z", " ", "\[Sigma]y"}]}], ")"}]}], "+", 
     RowBox[{"\[Alpha]", " ", "\[Beta]", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"Y", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
        RowBox[{"z", " ", 
         RowBox[{"(", 
          RowBox[{"\[Beta]", "+", 
           RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
       ")"}]}]}], 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"q", " ", "z"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"4", " ", "q", " ", "z"}], "-", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]]}], "}"}]], "Output",
 CellChangeTimes->{3.9623231443534794`*^9, 3.962323715597878*^9, 
  3.9623241277352257`*^9, 3.9625338313378944`*^9, 3.962597848791939*^9, 
  3.962600051320894*^9, 3.9713415115210533`*^9},
 CellLabel->"Out[38]=",ExpressionUUID->"cf9dd6cf-8edd-b141-8596-780e65f5b6a9"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\[Sigma]x,\[Sigma]y Lower Bound", "Section",
 CellChangeTimes->{{3.9623246574937305`*^9, 3.9623246996151657`*^9}, {
  3.962597773973393*^9, 
  3.962597786244307*^9}},ExpressionUUID->"0b4a5684-4873-3d44-bda0-\
f14e89ac48a7"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{
     RowBox[{"\[Sigma]x", "'"}], "[", 
     RowBox[{
     "\[Alpha]_", ",", "Y_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", 
      "L_"}], "]"}], ",", 
    RowBox[{
     RowBox[{"\[Sigma]y", "'"}], "[", 
     RowBox[{
     "\[Alpha]_", ",", "Y_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", 
      "L_"}], "]"}]}], "}"}], "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"\[Sigma]x", ",", "\[Sigma]y"}], "}"}], "/.", 
   RowBox[{"FullSimplify", "[", 
    RowBox[{
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"\[Sigma]x", ",", "\[Alpha]", ",", 
           RowBox[{
            RowBox[{
             RowBox[{"yAhat", "''"}], "[", 
             RowBox[{
             "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
              "z", ",", "q", ",", "\[Beta]", ",", "L"}], "]"}], "+", "0"}], 
           ",", 
           RowBox[{"pAeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", "z", ",", 
           RowBox[{
            RowBox[{"xAhat", "''"}], "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L"}], "]"}]}], "]"}], "==", 
         "0"}], "&&", 
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"\[Sigma]x", ",", "\[Alpha]", ",", 
           RowBox[{"Y", "-", 
            RowBox[{
             RowBox[{"yBhat", "''"}], "[", 
             RowBox[{
             "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
              "z", ",", "q", ",", "\[Beta]", ",", "L"}], "]"}], "+", "0"}], 
           ",", 
           RowBox[{"pBeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", "z", ",", 
           RowBox[{"1", "-", 
            RowBox[{
             RowBox[{"xBhat", "''"}], "[", 
             RowBox[{
             "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
              "z", ",", "q", ",", "\[Beta]", ",", "L"}], "]"}]}]}], "]"}], "==",
          "0"}], "&&", 
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"\[Sigma]y", ",", "\[Beta]", ",", 
           RowBox[{
            RowBox[{
             RowBox[{"xAhat", "''"}], "[", 
             RowBox[{
             "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
              "z", ",", "q", ",", "\[Beta]", ",", "L"}], "]"}], "+", "L"}], 
           ",", 
           RowBox[{"rAeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", "q", ",", 
           RowBox[{
            RowBox[{"yAhat", "''"}], "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L"}], "]"}]}], "]"}], "==", 
         "0"}], "&&", 
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"\[Sigma]y", ",", "\[Beta]", ",", 
           RowBox[{"1", "-", 
            RowBox[{
             RowBox[{"xBhat", "''"}], "[", 
             RowBox[{
             "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
              "z", ",", "q", ",", "\[Beta]", ",", "L"}], "]"}], "+", "0"}], 
           ",", 
           RowBox[{"rBeq", "[", 
            RowBox[{
            "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
             "z", ",", "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",",
              "0"}], "]"}], ",", "q", ",", 
           RowBox[{"Y", "-", 
            RowBox[{
             RowBox[{"yBhat", "''"}], "[", 
             RowBox[{
             "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", 
              "z", ",", "q", ",", "\[Beta]", ",", "L"}], "]"}]}]}], "]"}], "==",
          "0"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Sigma]x", ",", "\[Sigma]y"}], "}"}]}], "]"}], "[", 
     RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.962323176595495*^9, 3.962323547534401*^9}, {
  3.9623241332621403`*^9, 3.962324138177599*^9}, {3.962324172626326*^9, 
  3.96232420029942*^9}, {3.962324230798752*^9, 3.962324239889126*^9}, {
  3.962324421354349*^9, 3.9623244264734993`*^9}, {3.962600020682804*^9, 
  3.962600062370556*^9}},
 CellLabel->"In[39]:=",ExpressionUUID->"2ebbe932-4037-4f44-ba81-65c14d504b23"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", 
    FractionBox[
     RowBox[{
      RowBox[{"-", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"2", "+", "L"}], ")"}], " ", "q", " ", "z"}], ")"}]}], "+", 
      RowBox[{"q", " ", "Y", " ", "\[Alpha]"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "+", "L"}], ")"}], " ", "\[Alpha]", " ", "\[Beta]"}]}], 
     RowBox[{"2", " ", "q"}]]}], ",", 
   RowBox[{
    RowBox[{"q", " ", "Y"}], "-", 
    FractionBox[
     RowBox[{
      RowBox[{"(", 
       RowBox[{"z", "+", 
        RowBox[{"L", " ", "z"}], "+", 
        RowBox[{"Y", " ", "\[Alpha]"}]}], ")"}], " ", "\[Beta]"}], 
     RowBox[{"2", " ", "z"}]]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.962323522959276*^9, 3.9623235384621716`*^9}, 
   3.962323590049147*^9, 3.9623241439696693`*^9, 3.962324204021899*^9, 
   3.962324243176081*^9, 3.962324428669319*^9, 3.9625338317704296`*^9, 
   3.962597848947443*^9, 3.962600065556324*^9, 3.9713415116729317`*^9},
 CellLabel->"Out[39]=",ExpressionUUID->"1cbb047e-f000-684d-93f2-01bfecb71da4"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Utility Functions only include \
\[Sigma]x,\[Sigma]y,\[Alpha],\[Beta],z,q,Y,L", "Section",
 CellChangeTimes->{{3.962597813786089*^9, 3.962597835705454*^9}, 
   3.9625999202649555`*^9},ExpressionUUID->"40806356-fd71-5b41-a5d6-\
56bd84606337"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"u", "[", 
   RowBox[{"\[Sigma]x", ",", "\[Alpha]", ",", 
    RowBox[{
     RowBox[{
      RowBox[{"yAhat", "''"}], "[", 
      RowBox[{
      "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z", ",",
        "q", ",", "\[Beta]", ",", "L"}], "]"}], "+", "0"}], ",", 
    RowBox[{"pAeq", "[", 
     RowBox[{
     "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z", ",", 
      "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",", "0"}], "]"}], 
    ",", "z", ",", 
    RowBox[{
     RowBox[{"xAhat", "''"}], "[", 
     RowBox[{
     "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z", ",", 
      "q", ",", "\[Beta]", ",", "L"}], "]"}]}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"u", 
    RowBox[{"(", 
     RowBox[{"x", ";", "A"}], ")"}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"u", "[", 
   RowBox[{"\[Sigma]x", ",", "\[Alpha]", ",", 
    RowBox[{"Y", "-", 
     RowBox[{
      RowBox[{"yBhat", "''"}], "[", 
      RowBox[{
      "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z", ",",
        "q", ",", "\[Beta]", ",", "L"}], "]"}], "+", "0"}], ",", 
    RowBox[{"pBeq", "[", 
     RowBox[{
     "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z", ",", 
      "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",", "0"}], "]"}], 
    ",", "z", ",", 
    RowBox[{"1", "-", 
     RowBox[{
      RowBox[{"xBhat", "''"}], "[", 
      RowBox[{
      "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z", ",",
        "q", ",", "\[Beta]", ",", "L"}], "]"}]}]}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"u", 
    RowBox[{"(", 
     RowBox[{"x", ";", "B"}], ")"}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"u", "[", 
   RowBox[{"\[Sigma]y", ",", "\[Beta]", ",", 
    RowBox[{
     RowBox[{
      RowBox[{"xAhat", "''"}], "[", 
      RowBox[{
      "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z", ",",
        "q", ",", "\[Beta]", ",", "L"}], "]"}], "+", "L"}], ",", 
    RowBox[{"rAeq", "[", 
     RowBox[{
     "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z", ",", 
      "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",", "0"}], "]"}], 
    ",", "q", ",", 
    RowBox[{
     RowBox[{"yAhat", "''"}], "[", 
     RowBox[{
     "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z", ",", 
      "q", ",", "\[Beta]", ",", "L"}], "]"}]}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"u", 
    RowBox[{"(", 
     RowBox[{"y", ";", "A"}], ")"}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"u", "[", 
   RowBox[{"\[Sigma]y", ",", "\[Beta]", ",", 
    RowBox[{"1", "-", 
     RowBox[{
      RowBox[{"xBhat", "''"}], "[", 
      RowBox[{
      "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z", ",",
        "q", ",", "\[Beta]", ",", "L"}], "]"}], "+", "0"}], ",", 
    RowBox[{"rBeq", "[", 
     RowBox[{
     "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z", ",", 
      "q", ",", "\[Beta]", ",", "L", ",", "0", ",", "0", ",", "0"}], "]"}], 
    ",", "q", ",", 
    RowBox[{"Y", "-", 
     RowBox[{
      RowBox[{"yBhat", "''"}], "[", 
      RowBox[{
      "\[Sigma]x", ",", "\[Sigma]y", ",", "\[Alpha]", ",", "Y", ",", "z", ",",
        "q", ",", "\[Beta]", ",", "L"}], "]"}]}]}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"u", 
    RowBox[{"(", 
     RowBox[{"y", ";", "B"}], ")"}]}], 
   "*)"}]}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.962534189605047*^9, 3.96253425044363*^9}, {
   3.962534317645808*^9, 3.9625344035191593`*^9}, {3.962534454312462*^9, 
   3.96253449948666*^9}, 3.962534560276625*^9, {3.9626000685641346`*^9, 
   3.962600075900921*^9}},
 CellLabel->"In[40]:=",ExpressionUUID->"747f4495-10dc-6349-810f-acffd2fa190d"],

Cell[BoxData[
 RowBox[{"\[Sigma]x", "+", 
  FractionBox[
   RowBox[{
    RowBox[{
     RowBox[{"-", "2"}], " ", "q", " ", "Y", " ", "z", " ", "\[Alpha]"}], "+", 
    RowBox[{"z", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
    RowBox[{"Y", " ", 
     SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]"}], "+", 
    RowBox[{"2", " ", "L", " ", "z", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "q"}], " ", "z"}], "+", 
       RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "-", 
    RowBox[{"2", " ", "q", " ", "z", " ", "\[Sigma]x"}], "+", 
    RowBox[{"\[Alpha]", " ", "\[Beta]", " ", "\[Sigma]x"}], "+", 
    RowBox[{"z", " ", "\[Alpha]", " ", "\[Sigma]y"}]}], 
   RowBox[{
    RowBox[{"4", " ", "q", " ", "z"}], "-", 
    RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]], "+", 
  FractionBox[
   RowBox[{"z", " ", "\[Alpha]", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"-", "q"}], " ", "\[Beta]", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "+", 
      RowBox[{"2", " ", "q", " ", "z", " ", 
       RowBox[{"(", 
        RowBox[{"\[Beta]", "+", 
         RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}], "-", 
      RowBox[{"\[Alpha]", " ", "\[Beta]", " ", 
       RowBox[{"(", 
        RowBox[{"\[Beta]", "+", 
         RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
     ")"}]}], 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"q", " ", "z"}], "-", 
      RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", "q", " ", "z"}], "-", 
      RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]], "-", 
  FractionBox[
   RowBox[{"z", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"-", "L"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"q", " ", "z"}], "-", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", " ", "q", " ", "z"}], "-", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "+", 
      RowBox[{"q", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", " ", "q", " ", "z", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "-", 
         RowBox[{"\[Alpha]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"\[Beta]", " ", 
             RowBox[{"(", 
              RowBox[{"z", "+", 
               RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], 
            "+", 
            RowBox[{"z", " ", "\[Sigma]y"}]}], ")"}]}]}], ")"}]}]}], ")"}]}], 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"q", " ", "z"}], "-", 
      RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", "q", " ", "z"}], "-", 
      RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]]}]], "Output",
 CellChangeTimes->{3.9625345036258698`*^9, 3.962534579285124*^9, 
  3.962597849001873*^9, 3.9626000788655205`*^9, 3.9713415117205677`*^9},
 CellLabel->"Out[40]=",ExpressionUUID->"345ca96e-be52-af47-9452-f980d3d1d221"],

Cell[BoxData[
 RowBox[{"\[Sigma]x", "+", 
  FractionBox[
   RowBox[{
    RowBox[{
     RowBox[{"-", "2"}], " ", "q", " ", "Y", " ", "z", " ", "\[Alpha]"}], "+", 
    RowBox[{"z", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
    RowBox[{"L", " ", "z", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
    RowBox[{"Y", " ", 
     SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]"}], "-", 
    RowBox[{"2", " ", "q", " ", "z", " ", "\[Sigma]x"}], "+", 
    RowBox[{"\[Alpha]", " ", "\[Beta]", " ", "\[Sigma]x"}], "+", 
    RowBox[{"z", " ", "\[Alpha]", " ", "\[Sigma]y"}]}], 
   RowBox[{
    RowBox[{"4", " ", "q", " ", "z"}], "-", 
    RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]], "-", 
  RowBox[{"z", " ", 
   RowBox[{"(", 
    RowBox[{"1", "-", 
     FractionBox[
      RowBox[{
       RowBox[{
        SuperscriptBox["\[Alpha]", "2"], " ", 
        SuperscriptBox["\[Beta]", "2"]}], "+", 
       RowBox[{"2", " ", 
        SuperscriptBox["q", "2"], " ", "z", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"2", " ", "z"}], "-", 
          RowBox[{"Y", " ", "\[Alpha]"}], "-", "\[Sigma]x"}], ")"}]}], "+", 
       RowBox[{"q", " ", "\[Alpha]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"\[Beta]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "4"}], "+", "L"}], ")"}], " ", "z"}], "+", 
             RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "+", 
          RowBox[{"z", " ", "\[Sigma]y"}]}], ")"}]}]}], 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"q", " ", "z"}], "-", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", "q", " ", "z"}], "-", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]]}], ")"}]}], "+", 
  RowBox[{"\[Alpha]", " ", 
   RowBox[{"(", 
    RowBox[{"Y", "-", 
     FractionBox[
      RowBox[{
       RowBox[{"4", " ", 
        SuperscriptBox["q", "2"], " ", "Y", " ", 
        SuperscriptBox["z", "2"]}], "-", 
       RowBox[{"q", " ", "z", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"\[Beta]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{"2", "+", "L"}], ")"}], " ", "z"}], "+", 
             RowBox[{"4", " ", "Y", " ", "\[Alpha]"}], "-", "\[Sigma]x"}], 
            ")"}]}], "+", 
          RowBox[{"2", " ", "z", " ", "\[Sigma]y"}]}], ")"}]}], "+", 
       RowBox[{"\[Alpha]", " ", "\[Beta]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Y", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
          RowBox[{"z", " ", 
           RowBox[{"(", 
            RowBox[{"\[Beta]", "+", 
             RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
         ")"}]}]}], 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"q", " ", "z"}], "-", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", "q", " ", "z"}], "-", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]]}], 
    ")"}]}]}]], "Output",
 CellChangeTimes->{3.9625345036258698`*^9, 3.962534579285124*^9, 
  3.962597849001873*^9, 3.9626000788655205`*^9, 3.9713415117205677`*^9},
 CellLabel->"Out[41]=",ExpressionUUID->"7e48b32b-1706-2347-816e-b144472b921b"],

Cell[BoxData[
 RowBox[{"\[Sigma]y", "-", 
  FractionBox[
   RowBox[{"q", " ", "z", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"-", "q"}], " ", "\[Beta]", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "+", 
      RowBox[{"2", " ", "q", " ", "z", " ", 
       RowBox[{"(", 
        RowBox[{"\[Beta]", "+", 
         RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}], "-", 
      RowBox[{"\[Alpha]", " ", "\[Beta]", " ", 
       RowBox[{"(", 
        RowBox[{"\[Beta]", "+", 
         RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
     ")"}]}], 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"q", " ", "z"}], "-", 
      RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", "q", " ", "z"}], "-", 
      RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]], "-", 
  FractionBox[
   RowBox[{
    RowBox[{
     RowBox[{"-", "\[Alpha]"}], " ", "\[Beta]", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "+", "L"}], ")"}], " ", "\[Beta]"}], "+", "\[Sigma]y"}],
       ")"}]}], "+", 
    RowBox[{"q", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"2", "+", 
          RowBox[{"2", " ", "L"}]}], ")"}], " ", "z", " ", "\[Beta]"}], "-", 
       RowBox[{"\[Beta]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "+", 
       RowBox[{"2", " ", "z", " ", "\[Sigma]y"}]}], ")"}]}]}], 
   RowBox[{
    RowBox[{"4", " ", "q", " ", "z"}], "-", 
    RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]], "+", 
  RowBox[{"\[Beta]", " ", 
   RowBox[{"(", 
    RowBox[{"L", "+", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", "L"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"q", " ", "z"}], "-", 
          RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"2", " ", "q", " ", "z"}], "-", 
          RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}], "+", 
       RowBox[{"q", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"2", " ", "q", " ", "z", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "-", 
          RowBox[{"\[Alpha]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"\[Beta]", " ", 
              RowBox[{"(", 
               RowBox[{"z", "+", 
                RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], 
             "+", 
             RowBox[{"z", " ", "\[Sigma]y"}]}], ")"}]}]}], ")"}]}]}], 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"q", " ", "z"}], "-", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", "q", " ", "z"}], "-", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]]}], 
    ")"}]}]}]], "Output",
 CellChangeTimes->{3.9625345036258698`*^9, 3.962534579285124*^9, 
  3.962597849001873*^9, 3.9626000788655205`*^9, 3.97134151173497*^9},
 CellLabel->"Out[42]=",ExpressionUUID->"3a0560b4-1d04-8e4d-9555-fa7cb7aca224"],

Cell[BoxData[
 RowBox[{"\[Sigma]y", "-", 
  FractionBox[
   RowBox[{
    RowBox[{
     RowBox[{"-", "\[Alpha]"}], " ", "\[Beta]", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "+", "L"}], ")"}], " ", "\[Beta]"}], "+", "\[Sigma]y"}],
       ")"}]}], "+", 
    RowBox[{"q", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"2", "+", "L"}], ")"}], " ", "z", " ", "\[Beta]"}], "-", 
       RowBox[{"\[Beta]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "+", 
       RowBox[{"2", " ", "z", " ", "\[Sigma]y"}]}], ")"}]}]}], 
   RowBox[{
    RowBox[{"4", " ", "q", " ", "z"}], "-", 
    RowBox[{"\[Alpha]", " ", "\[Beta]"}]}]], "+", 
  RowBox[{"\[Beta]", " ", 
   RowBox[{"(", 
    RowBox[{"1", "-", 
     FractionBox[
      RowBox[{
       RowBox[{
        SuperscriptBox["\[Alpha]", "2"], " ", 
        SuperscriptBox["\[Beta]", "2"]}], "+", 
       RowBox[{"2", " ", 
        SuperscriptBox["q", "2"], " ", "z", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"2", " ", "z"}], "-", 
          RowBox[{"Y", " ", "\[Alpha]"}], "-", "\[Sigma]x"}], ")"}]}], "+", 
       RowBox[{"q", " ", "\[Alpha]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"\[Beta]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "4"}], "+", "L"}], ")"}], " ", "z"}], "+", 
             RowBox[{"Y", " ", "\[Alpha]"}], "+", "\[Sigma]x"}], ")"}]}], "+", 
          RowBox[{"z", " ", "\[Sigma]y"}]}], ")"}]}]}], 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"q", " ", "z"}], "-", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", "q", " ", "z"}], "-", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]]}], ")"}]}], "-", 
  RowBox[{"q", " ", 
   RowBox[{"(", 
    RowBox[{"Y", "-", 
     FractionBox[
      RowBox[{
       RowBox[{"4", " ", 
        SuperscriptBox["q", "2"], " ", "Y", " ", 
        SuperscriptBox["z", "2"]}], "-", 
       RowBox[{"q", " ", "z", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"\[Beta]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{"2", "+", "L"}], ")"}], " ", "z"}], "+", 
             RowBox[{"4", " ", "Y", " ", "\[Alpha]"}], "-", "\[Sigma]x"}], 
            ")"}]}], "+", 
          RowBox[{"2", " ", "z", " ", "\[Sigma]y"}]}], ")"}]}], "+", 
       RowBox[{"\[Alpha]", " ", "\[Beta]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Y", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
          RowBox[{"z", " ", 
           RowBox[{"(", 
            RowBox[{"\[Beta]", "+", 
             RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}]}], 
         ")"}]}]}], 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"q", " ", "z"}], "-", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", "q", " ", "z"}], "-", 
         RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}]}]]}], 
    ")"}]}]}]], "Output",
 CellChangeTimes->{3.9625345036258698`*^9, 3.962534579285124*^9, 
  3.962597849001873*^9, 3.9626000788655205`*^9, 3.97134151174115*^9},
 CellLabel->"Out[43]=",ExpressionUUID->"f2bc6e6b-c162-f446-88f2-297252e4cf48"]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Find Subsidies", "Chapter",
 CellChangeTimes->{{3.9625986905493317`*^9, 
  3.962598701809437*^9}},ExpressionUUID->"4945309f-a7de-5f42-9796-\
cd1dd11225cd"],

Cell[CellGroupData[{

Cell["Specific Value to test Subsidies", "Section",
 CellChangeTimes->{{3.962529117367857*^9, 3.962529122831496*^9}, {
  3.9625980253813686`*^9, 
  3.962598051122851*^9}},ExpressionUUID->"abf52c42-3e69-8b41-bb64-\
f9463486e6d8"],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Sigma]x1", "=", "0.4"}], ";", 
  RowBox[{"\[Sigma]y1", "=", "1.5"}], ";", 
  RowBox[{"\[Alpha]1", "=", "0.7"}], ";", 
  RowBox[{"Y1", "=", "1.6"}], ";", 
  RowBox[{"z1", "=", "0.9"}], ";", 
  RowBox[{"q1", "=", "1.4"}], ";", 
  RowBox[{"\[Beta]1", "=", "1.1"}], ";", 
  RowBox[{"LxA1", "=", "1.6"}], ";", 
  RowBox[{"LxB", "=", "0"}], ";", 
  RowBox[{"LyA", "=", "0"}], ";", 
  RowBox[{"LyB", "=", "0"}], ";"}]], "Input",
 CellChangeTimes->{{3.9624018547494392`*^9, 3.96240186729871*^9}, 
   3.962402818779262*^9, {3.962439678786089*^9, 3.9624397125700912`*^9}, {
   3.9624409438016567`*^9, 3.962440956326483*^9}, {3.962441086270864*^9, 
   3.962441086399193*^9}, {3.962441135710407*^9, 3.9624411616251564`*^9}, {
   3.962441210720781*^9, 3.9624412389865932`*^9}, {3.9624412864676285`*^9, 
   3.962441342641516*^9}, 3.962441402279466*^9, 3.962441439370852*^9, 
   3.962441502101073*^9, 3.962441553949335*^9, 3.962441595263168*^9, 
   3.962441628179287*^9, {3.9624416704878063`*^9, 3.96244167977742*^9}, {
   3.962441723902429*^9, 3.962441742405972*^9}, 3.962441774150894*^9, {
   3.962441816146927*^9, 3.962441819939178*^9}, {3.9624418571091747`*^9, 
   3.962441938619362*^9}, {3.962442128022737*^9, 3.962442177064535*^9}, 
   3.962442208685318*^9, {3.962442279351661*^9, 3.962442303663152*^9}, {
   3.962442351769159*^9, 3.962442402286068*^9}, 3.962442432386102*^9, {
   3.962442563263733*^9, 3.9624425688407135`*^9}, 3.962442643119838*^9, 
   3.9624427338824997`*^9, 3.9624428610227566`*^9, 3.962442958030342*^9, 
   3.962443172560799*^9, 3.962443249603508*^9, 3.962443419248834*^9, 
   3.962443516463449*^9, 3.962443610873655*^9, 3.962443709318771*^9, {
   3.962443918772955*^9, 3.962443919322882*^9}, {3.962443974362888*^9, 
   3.962444000483677*^9}, 3.9624440599909*^9, {3.962528205131932*^9, 
   3.962528208218586*^9}, {3.962528258164055*^9, 3.962528259712384*^9}, 
   3.962528303009552*^9, {3.962528336855339*^9, 3.962528361367201*^9}, 
   3.962528399276993*^9, 3.962528436038894*^9, {3.9625284724533863`*^9, 
   3.9625285012402935`*^9}, 3.96252893702326*^9, {3.962528973183468*^9, 
   3.96252898365468*^9}, {3.962530495025591*^9, 3.962530508750183*^9}, {
   3.962530549467009*^9, 3.96253057863546*^9}, 3.96253075023283*^9, {
   3.9625309412388787`*^9, 3.9625309606235065`*^9}, 3.962531028516535*^9, 
   3.9625310859290295`*^9, 3.962531146960238*^9, {3.962531193434635*^9, 
   3.9625312181957912`*^9}, 3.9625312892898254`*^9, 3.9625313897883205`*^9, 
   3.9625314239266148`*^9, {3.962531477345188*^9, 3.962531498044327*^9}, {
   3.962531531824221*^9, 3.9625315545749016`*^9}, {3.9625316200476875`*^9, 
   3.96253162021167*^9}, {3.962531681722138*^9, 3.96253175479941*^9}, {
   3.962531791952832*^9, 3.962531823823273*^9}, {3.962531855810898*^9, 
   3.962531871067965*^9}, 3.962532093239092*^9, {3.962532127255167*^9, 
   3.96253229890629*^9}, {3.962532335273117*^9, 3.962532362035258*^9}, {
   3.9625324021072483`*^9, 3.962532426780306*^9}, {3.962532548820612*^9, 
   3.962532557509627*^9}, 3.962532588836066*^9, 3.9625326615494995`*^9, {
   3.962532718444758*^9, 3.9625327259808826`*^9}, {3.962532843156868*^9, 
   3.9625329265235043`*^9}, {3.962532968196383*^9, 3.962532989579643*^9}},
 CellLabel->"In[44]:=",ExpressionUUID->"5fb78389-c9f4-6c45-914d-179ce1329433"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
    RowBox[{"pA", "'"}], ",", 
    RowBox[{"rA", "'"}], ",", 
    RowBox[{"pB", "'"}], ",", 
    RowBox[{"rB", "'"}]}], "}"}], "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"pA", ",", "rA", ",", "pB", ",", "rB"}], "}"}], "/.", 
   RowBox[{
    RowBox[{"FullSimplify", "[", 
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"pA", "==", 
         RowBox[{"pAeq", "[", 
          RowBox[{
          "\[Sigma]x1", ",", "\[Sigma]y1", ",", "\[Alpha]1", ",", "Y1", ",", 
           "z1", ",", "q1", ",", "\[Beta]1", ",", "LxA1", ",", "LxB", ",", 
           "LyA", ",", "LyB"}], "]"}]}], "&&", 
        RowBox[{"rA", "==", 
         RowBox[{"rAeq", "[", 
          RowBox[{
          "\[Sigma]x1", ",", "\[Sigma]y1", ",", "\[Alpha]1", ",", "Y1", ",", 
           "z1", ",", "q1", ",", "\[Beta]1", ",", "LxA1", ",", "LxB", ",", 
           "LyA", ",", "LyB"}], "]"}]}], "&&", 
        RowBox[{"pB", "==", 
         RowBox[{"pBeq", "[", 
          RowBox[{
          "\[Sigma]x1", ",", "\[Sigma]y1", ",", "\[Alpha]1", ",", "Y1", ",", 
           "z1", ",", "q1", ",", "\[Beta]1", ",", "LxA1", ",", "LxB", ",", 
           "LyA", ",", "LyB"}], "]"}]}], "&&", 
        RowBox[{"rB", "==", 
         RowBox[{"rBeq", "[", 
          RowBox[{
          "\[Sigma]x1", ",", "\[Sigma]y1", ",", "\[Alpha]1", ",", "Y1", ",", 
           "z1", ",", "q1", ",", "\[Beta]1", ",", "LxA1", ",", "LxB", ",", 
           "LyA", ",", "LyB"}], "]"}]}]}], ",", 
       RowBox[{"{", 
        RowBox[{"pA", ",", "rA", ",", "pB", ",", "rB"}], "}"}]}], "]"}], 
     "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.962402703253011*^9, 3.962402726581137*^9}, {
  3.9624027851099434`*^9, 3.962402864999265*^9}, {3.962439726670042*^9, 
  3.962439755188753*^9}},
 CellLabel->"In[45]:=",ExpressionUUID->"90ff39bb-5d09-c644-b727-ea98751d0791"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.5698360655737704`", ",", "1.2386885245901642`", ",", 
   RowBox[{"-", "0.02032786885245917`"}], ",", "0.719344262295082`"}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.962529267142611*^9, 3.962529278839722*^9}, {
   3.962530503886566*^9, 3.9625305129430084`*^9}, 3.962530583287974*^9, 
   3.962531089962916*^9, 3.962531150374012*^9, 3.962531368122614*^9, {
   3.9625314030327225`*^9, 3.962531428047177*^9}, {3.962531481562929*^9, 
   3.96253150149523*^9}, {3.9625315349351044`*^9, 3.9625315574592*^9}, 
   3.962531623501732*^9, {3.962531685284029*^9, 3.962531713410591*^9}, 
   3.9625317581739*^9, {3.962531861580755*^9, 3.9625318742948685`*^9}, 
   3.962532097326456*^9, {3.962532131545749*^9, 3.9625321551522617`*^9}, {
   3.962532187449383*^9, 3.962532279336258*^9}, {3.962532339563305*^9, 
   3.9625323660339775`*^9}, 3.962532406853358*^9, {3.962532545760088*^9, 
   3.962532561374344*^9}, 3.962532676628332*^9, {3.962532722516604*^9, 
   3.962532730797175*^9}, {3.962532847654207*^9, 3.962532872977274*^9}, {
   3.962532923954305*^9, 3.962532930532076*^9}, 3.962533010287977*^9, 
   3.962598115825687*^9, 3.962598297727991*^9, 3.9625995892286377`*^9, 
   3.9626000885340233`*^9, 3.9713415188858128`*^9},
 CellLabel->"Out[45]=",ExpressionUUID->"6e9bd15f-2e4b-8841-9c12-d07ed6d5c49a"]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"f\[Sigma]x", "[", 
    RowBox[{
    "\[Alpha]_", ",", "Y_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", 
     "L_"}], "]"}], ":=", 
   RowBox[{"-", 
    FractionBox[
     RowBox[{
      RowBox[{"-", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"2", "+", "L"}], ")"}], " ", "q", " ", "z"}], ")"}]}], "+", 
      RowBox[{"q", " ", "Y", " ", "\[Alpha]"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "+", "L"}], ")"}], " ", "\[Alpha]", " ", "\[Beta]"}]}], 
     RowBox[{"2", " ", "q"}]]}]}], 
  RowBox[{"(*", 
   RowBox[{"\[Sigma]x", ">=", 
    RowBox[{"-", 
     FractionBox[
      RowBox[{
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"2", "+", "L"}], ")"}], " ", "q", " ", "z"}], ")"}]}], "+", 
       RowBox[{"q", " ", "Y", " ", "\[Alpha]"}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "+", "L"}], ")"}], " ", "\[Alpha]", " ", "\[Beta]"}]}], 
      RowBox[{"2", " ", "q"}]]}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"f\[Sigma]y", "[", 
   RowBox[{
   "\[Alpha]_", ",", "Y_", ",", "z_", ",", "q_", ",", "\[Beta]_", ",", "L_"}],
    "]"}], ":=", 
  RowBox[{
   RowBox[{"q", " ", "Y"}], "-", 
   FractionBox[
    RowBox[{
     RowBox[{"(", 
      RowBox[{"z", "+", 
       RowBox[{"L", " ", "z"}], "+", 
       RowBox[{"Y", " ", "\[Alpha]"}]}], ")"}], " ", "\[Beta]"}], 
    RowBox[{"2", " ", "z"}]], 
   RowBox[{"(*", 
    RowBox[{"\[Sigma]y", ">=", 
     RowBox[{
      RowBox[{"q", " ", "Y"}], "-", 
      FractionBox[
       RowBox[{
        RowBox[{"(", 
         RowBox[{"z", "+", 
          RowBox[{"L", " ", "z"}], "+", 
          RowBox[{"Y", " ", "\[Alpha]"}]}], ")"}], " ", "\[Beta]"}], 
       RowBox[{"2", " ", "z"}]]}]}], "*)"}]}]}]}], "Input",
 CellChangeTimes->{{3.96244003390854*^9, 3.962440233442196*^9}, {
  3.9624402686366653`*^9, 3.962440277775135*^9}, {3.9624403208173738`*^9, 
  3.962440323525533*^9}, {3.9624404011499367`*^9, 3.962440412139705*^9}, {
  3.962440526207878*^9, 3.962440561601555*^9}},
 CellLabel->"In[46]:=",ExpressionUUID->"002a7799-b916-2249-ab40-b7777a612d04"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"f\[Sigma]x", "[", 
   RowBox[{
   "\[Alpha]1", ",", "Y1", ",", "z1", ",", "q1", ",", "\[Beta]1", ",", 
    "LxA1"}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"\[Sigma]x", " ", "Lower", " ", "Bound"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.962440428184698*^9, 3.962440482214124*^9}, {
  3.962529180562744*^9, 3.962529197236479*^9}, {3.9625980707171288`*^9, 
  3.9625980818100834`*^9}},
 CellLabel->"In[48]:=",ExpressionUUID->"9141a280-bbf2-7040-8c29-1a18e9f8c7bb"],

Cell[BoxData["0.3450000000000004`"], "Output",
 CellChangeTimes->{
  3.962440483075283*^9, 3.9624405728639736`*^9, 3.962440972015205*^9, 
   3.962441097979599*^9, 3.962441172588892*^9, 3.96244124688732*^9, 
   3.9624413664523506`*^9, 3.962441410290066*^9, 3.962441525110897*^9, 
   3.962441603354818*^9, 3.9624416361325474`*^9, 3.9624416884785423`*^9, {
   3.9624417327039585`*^9, 3.962441750318531*^9}, 3.9624418714735928`*^9, 
   3.962441956298506*^9, 3.962442136687151*^9, 3.9624423607917767`*^9, {
   3.962442395500944*^9, 3.962442409129723*^9}, 3.962442705753648*^9, 
   3.9624428105817347`*^9, 3.962442869523186*^9, 3.962443062740282*^9, 
   3.962443329816627*^9, 3.962443484391897*^9, 3.962443576193964*^9, 
   3.962443664003828*^9, 3.962444074776109*^9, 3.962528231494484*^9, 
   3.962528311516466*^9, 3.9625283454157314`*^9, {3.962528380543197*^9, 
   3.9625284087401543`*^9}, 3.9625284433417664`*^9, {3.962528480288109*^9, 
   3.962528509354854*^9}, 3.9625289522556705`*^9, 3.9625289897895966`*^9, 
   3.962529285356327*^9, 3.962530521145748*^9, 3.962530587728459*^9, 
   3.9625310952836876`*^9, 3.962531155141802*^9, {3.9625314133915653`*^9, 
   3.962531430931078*^9}, {3.962531483905407*^9, 3.962531503982624*^9}, {
   3.962531537199892*^9, 3.9625315595284805`*^9}, 3.962531627127054*^9, {
   3.96253169672106*^9, 3.962531715826227*^9}, 3.962531771576475*^9, 
   3.962532100630806*^9, {3.962532134701177*^9, 3.962532158236334*^9}, {
   3.962532223379881*^9, 3.96253224505564*^9}, 3.962532282608879*^9, {
   3.962532348023924*^9, 3.9625323690154057`*^9}, 3.9625324121200447`*^9, 
   3.962532563643473*^9, 3.9625326816721687`*^9, 3.962532936340217*^9, 
   3.962533012798132*^9, 3.96259811585639*^9, 3.962598297755993*^9, 
   3.96259958926457*^9, 3.962600088568651*^9, 3.9713415189178677`*^9},
 CellLabel->"Out[48]=",ExpressionUUID->"c02a5c13-1f18-5346-820e-8dedf1e79bf1"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"f\[Sigma]y", "[", 
   RowBox[{
   "\[Alpha]1", ",", "Y1", ",", "z1", ",", "q1", ",", "\[Beta]1", ",", 
    "LxA1"}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"\[Sigma]y", " ", "Lower", " ", "Bound"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.962440491194668*^9, 3.9624405005667553`*^9}, {
  3.962529199437687*^9, 3.9625292209728355`*^9}, {3.962598084723905*^9, 
  3.9625980900566597`*^9}},
 CellLabel->"In[49]:=",ExpressionUUID->"244b78b3-f103-da4a-85ff-469d432fcbf2"],

Cell[BoxData["0.1255555555555552`"], "Output",
 CellChangeTimes->{
  3.962440502193077*^9, 3.962440574932726*^9, 3.9624409762806168`*^9, 
   3.962441099545051*^9, 3.962441174203528*^9, 3.96244124934256*^9, 
   3.9624413695766983`*^9, 3.96244141274723*^9, 3.962441526398924*^9, 
   3.962441605434946*^9, 3.962441638140899*^9, 3.962441690392509*^9, {
   3.96244173459144*^9, 3.96244175219388*^9}, 3.962441872869478*^9, 
   3.962441957763187*^9, 3.962442138679958*^9, 3.962442184637424*^9, 
   3.962442362362247*^9, {3.962442396966531*^9, 3.962442410779873*^9}, 
   3.962442707416881*^9, 3.962442819452713*^9, 3.9624428713349*^9, 
   3.9624430658103905`*^9, 3.96244338016612*^9, 3.96244348636959*^9, 
   3.9624435826131363`*^9, 3.962443676128093*^9, 3.9625282336726475`*^9, 
   3.9625283135318108`*^9, 3.9625283476112995`*^9, {3.962528382763443*^9, 
   3.962528411071945*^9}, 3.962528445494137*^9, 3.962528482652214*^9, 
   3.9625285181684284`*^9, 3.962528954176529*^9, 3.962528991729625*^9, 
   3.962529288179631*^9, 3.962530523962057*^9, 3.962530589872422*^9, 
   3.962531097044623*^9, 3.962531156854616*^9, {3.962531415675619*^9, 
   3.962531432927183*^9}, {3.962531485560011*^9, 3.962531505883116*^9}, {
   3.962531538888953*^9, 3.962531561221294*^9}, 3.962531629114014*^9, {
   3.962531698330483*^9, 3.962531717638348*^9}, 3.962531774392861*^9, 
   3.962532102725489*^9, {3.962532136666509*^9, 3.9625321603708344`*^9}, {
   3.962532225424871*^9, 3.962532247280367*^9}, 3.962532284798067*^9, {
   3.9625323500818653`*^9, 3.962532371598072*^9}, 3.962532414531271*^9, 
   3.9625325655539417`*^9, 3.96253268346039*^9, 3.962532938189968*^9, 
   3.9625330146631203`*^9, 3.9625981158828278`*^9, 3.962598297786066*^9, 
   3.962599589295233*^9, 3.962600088598823*^9, 3.9713415189518604`*^9},
 CellLabel->"Out[49]=",ExpressionUUID->"de9d5bc4-08fe-1344-b03c-6161d2a9d247"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"qz", ">", "\[Alpha]\[Beta]"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"fqz", "[", 
    RowBox[{"\[Alpha]_", ",", "z_", ",", "q_", ",", "\[Beta]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"q", "*", "z"}], "-", 
    RowBox[{"\[Alpha]", "*", "\[Beta]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.962440594599079*^9, 3.9624407122356434`*^9}},
 CellLabel->"In[50]:=",ExpressionUUID->"d8f4e21a-0f1e-1648-8a0a-e09550d07089"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"fqz", "[", 
   RowBox[{"\[Alpha]1", ",", "z1", ",", "q1", ",", "\[Beta]1"}], "]"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"qz", "-", "\[Alpha]\[Beta]"}], ">", "0"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.962440722877529*^9, 3.9624407258199863`*^9}, {
  3.962529300729847*^9, 3.962529322571324*^9}},
 CellLabel->"In[51]:=",ExpressionUUID->"90caf400-2533-b940-9809-a224eeef115b"],

Cell[BoxData["0.49`"], "Output",
 CellChangeTimes->{
  3.962440729187271*^9, 3.962440984812349*^9, 3.962441109918901*^9, 
   3.962441179025469*^9, 3.962441254455948*^9, 3.962441373769348*^9, 
   3.962441417277151*^9, 3.962441534086708*^9, 3.962441609734716*^9, 
   3.962441641399166*^9, 3.962441693742243*^9, 3.9624417571710453`*^9, 
   3.962441960979795*^9, 3.9624421867786*^9, 3.962442413283376*^9, 
   3.962442710603615*^9, 3.962442874169697*^9, 3.962443075746589*^9, 
   3.962443388661981*^9, 3.962443494071216*^9, 3.9624435878049183`*^9, 
   3.96244367797394*^9, 3.9625282391976223`*^9, 3.9625283171760464`*^9, 
   3.9625283506783257`*^9, 3.962528385592798*^9, 3.9625284168721905`*^9, 
   3.962528448630051*^9, 3.962528486231544*^9, {3.9625285241460037`*^9, 
   3.962528535755516*^9}, 3.962528959786289*^9, 3.962528993743955*^9, 
   3.962529305262057*^9, 3.962530528628247*^9, 3.962530593533174*^9, 
   3.962531099278923*^9, 3.9625311593792915`*^9, {3.962531417956303*^9, 
   3.962531436384508*^9}, {3.9625314876702595`*^9, 3.962531508370425*^9}, {
   3.962531541527842*^9, 3.962531563216507*^9}, 3.962531631081311*^9, {
   3.962531700628298*^9, 3.962531719468973*^9}, 3.962531776780875*^9, 
   3.9625321054295*^9, {3.962532139144943*^9, 3.962532162861326*^9}, {
   3.962532227454645*^9, 3.962532249290114*^9}, 3.962532287336024*^9, {
   3.962532352910368*^9, 3.962532373905855*^9}, 3.962532417314804*^9, 
   3.962532567974118*^9, 3.962532690049684*^9, 3.96253294151648*^9, 
   3.962533017131558*^9, 3.962598115916761*^9, 3.962598297816935*^9, 
   3.962599589332394*^9, 3.962600088634863*^9, 3.9713415189651146`*^9},
 CellLabel->"Out[51]=",ExpressionUUID->"f8d7c523-ac41-9f4e-a6a2-eefc5d4e81a7"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"\[Sigma]x", "<=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"64", " ", 
        SuperscriptBox["q", "4"], " ", 
        SuperscriptBox["z", "4"], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"L", " ", "z"}], "+", 
          RowBox[{"Y", " ", "\[Alpha]"}]}], ")"}]}], "+", 
       RowBox[{"8", " ", 
        SuperscriptBox["q", "2"], " ", 
        SuperscriptBox["z", "2"], " ", 
        SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"4", "+", 
             RowBox[{"11", " ", "L"}]}], ")"}], " ", "z", " ", "\[Beta]"}], 
          "+", 
          RowBox[{"7", " ", "Y", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
          RowBox[{"4", " ", "z", " ", "\[Sigma]y"}]}], ")"}]}], "+", 
       RowBox[{
        SuperscriptBox["\[Alpha]", "3"], " ", 
        SuperscriptBox["\[Beta]", "3"], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"\[Beta]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"2", " ", "\[Alpha]", " ", 
              RowBox[{"(", 
               RowBox[{"z", "+", 
                RowBox[{"2", " ", "L", " ", "z"}], "+", 
                RowBox[{"Y", " ", "\[Alpha]"}]}], ")"}]}], "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"1", "+", "L"}], ")"}], " ", "z", " ", "\[Beta]"}]}], 
            ")"}]}], "+", 
          RowBox[{"z", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"2", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}], " ", 
           "\[Sigma]y"}]}], ")"}]}], "-", 
       RowBox[{"32", " ", 
        SuperscriptBox["q", "3"], " ", 
        SuperscriptBox["z", "3"], " ", "\[Alpha]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"3", " ", "Y", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
          RowBox[{"z", " ", 
           RowBox[{"(", 
            RowBox[{"\[Beta]", "+", 
             RowBox[{"4", " ", "L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], 
            ")"}]}]}], ")"}]}], "-", 
       RowBox[{"q", " ", "z", " ", 
        SuperscriptBox["\[Alpha]", "2"], " ", 
        SuperscriptBox["\[Beta]", "2"], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Y", " ", "\[Alpha]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"16", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], " ", 
           "\[Beta]"}], "+", 
          RowBox[{"2", " ", "z", " ", "\[Beta]", " ", 
           RowBox[{"(", 
            RowBox[{"\[Beta]", "+", 
             RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}], "+", 
          RowBox[{"4", " ", "z", " ", "\[Alpha]", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{"3", "+", 
                RowBox[{"7", " ", "L"}]}], ")"}], " ", "\[Beta]"}], "+", 
             RowBox[{"3", " ", "\[Sigma]y"}]}], ")"}]}]}], ")"}]}], "-", 
       RowBox[{"2", " ", "q", " ", "z", " ", 
        RowBox[{"\[Sqrt]", 
         RowBox[{"(", 
          RowBox[{
           FractionBox["1", 
            RowBox[{
             SuperscriptBox["q", "2"], " ", 
             SuperscriptBox["z", "2"]}]], 
           SuperscriptBox["\[Alpha]", "2"], " ", 
           SuperscriptBox["\[Beta]", "2"], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"2", " ", "q", " ", "z"}], "-", 
             RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"2", " ", 
              SuperscriptBox["L", "2"], " ", 
              SuperscriptBox["z", "2"], " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"64", " ", 
                 SuperscriptBox["q", "5"], " ", 
                 SuperscriptBox["z", "5"]}], "+", 
                RowBox[{"4", " ", 
                 SuperscriptBox["q", "3"], " ", 
                 SuperscriptBox["z", "3"], " ", "\[Alpha]", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"53", " ", "\[Alpha]"}], "-", 
                   RowBox[{"2", " ", "\[Beta]"}]}], ")"}], " ", 
                 SuperscriptBox["\[Beta]", "2"]}], "+", 
                RowBox[{"q", " ", "z", " ", 
                 SuperscriptBox["\[Alpha]", "3"], " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"24", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                 " ", 
                 SuperscriptBox["\[Beta]", "4"]}], "-", 
                RowBox[{"2", " ", 
                 SuperscriptBox["\[Alpha]", "5"], " ", 
                 SuperscriptBox["\[Beta]", "5"]}], "+", 
                RowBox[{"4", " ", 
                 SuperscriptBox["q", "4"], " ", 
                 SuperscriptBox["z", "4"], " ", "\[Beta]", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"-", "48"}], " ", "\[Alpha]"}], "+", "\[Beta]"}], 
                  ")"}]}], "+", 
                RowBox[{"2", " ", 
                 SuperscriptBox["q", "2"], " ", 
                 SuperscriptBox["z", "2"], " ", 
                 SuperscriptBox["\[Alpha]", "2"], " ", 
                 SuperscriptBox["\[Beta]", "3"], " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"-", "53"}], " ", "\[Alpha]"}], "+", 
                   RowBox[{"3", " ", "\[Beta]"}]}], ")"}]}]}], ")"}]}], "+", 
             RowBox[{"L", " ", "z", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"4", " ", "q", " ", "z"}], "-", 
                RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"64", " ", 
                 SuperscriptBox["q", "4"], " ", "Y", " ", 
                 SuperscriptBox["z", "4"], " ", "\[Alpha]"}], "+", 
                RowBox[{"2", " ", 
                 SuperscriptBox["q", "2"], " ", 
                 SuperscriptBox["z", "2"], " ", "\[Alpha]", " ", "\[Beta]", 
                 " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"Y", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"56", " ", "\[Alpha]"}], "+", 
                    RowBox[{"5", " ", "\[Beta]"}]}], ")"}]}], "+", 
                   RowBox[{"2", " ", "z", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"14", " ", "\[Alpha]"}], "-", 
                    RowBox[{"3", " ", "\[Beta]"}]}], ")"}], " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], ")"}]}],
                 "-", 
                RowBox[{"2", " ", "q", " ", "z", " ", 
                 SuperscriptBox["\[Alpha]", "2"], " ", 
                 SuperscriptBox["\[Beta]", "2"], " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                   "3", " ", "Y", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"6", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}]}],
                    "+", 
                   RowBox[{"2", " ", "z", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"7", " ", "\[Alpha]"}], "-", 
                    RowBox[{"2", " ", "\[Beta]"}]}], ")"}], " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], ")"}]}],
                 "+", 
                RowBox[{
                 SuperscriptBox["\[Alpha]", "3"], " ", 
                 SuperscriptBox["\[Beta]", "3"], " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"Y", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"4", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}]}],
                    "+", 
                   RowBox[{"z", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"4", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                    " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], ")"}]}],
                 "-", 
                RowBox[{"8", " ", 
                 SuperscriptBox["q", "3"], " ", 
                 SuperscriptBox["z", "3"], " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"Y", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"18", " ", "\[Alpha]"}], "+", "\[Beta]"}], 
                    ")"}]}], "+", 
                   RowBox[{"z", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"4", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                    " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], 
                  ")"}]}]}], ")"}]}], "+", 
             RowBox[{
              SuperscriptBox[
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "4"}], " ", "q", " ", "z"}], "+", 
                 RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], "2"], " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"8", " ", 
                 SuperscriptBox["q", "3"], " ", 
                 SuperscriptBox["Y", "2"], " ", 
                 SuperscriptBox["z", "3"], " ", 
                 SuperscriptBox["\[Alpha]", "2"]}], "+", 
                RowBox[{"2", " ", 
                 SuperscriptBox["q", "2"], " ", 
                 SuperscriptBox["z", "2"], " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    SuperscriptBox["Y", "2"], " ", 
                    SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"-", "6"}], " ", "\[Alpha]"}], "+", "\[Beta]"}], 
                    ")"}]}], "-", 
                   RowBox[{"2", " ", "Y", " ", "z", " ", "\[Alpha]", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"2", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}], 
                    " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "+", 
                   RowBox[{
                    SuperscriptBox["z", "2"], " ", 
                    SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}]}], 
                  ")"}]}], "-", 
                RowBox[{
                 SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    SuperscriptBox["Y", "2"], " ", 
                    SuperscriptBox["\[Alpha]", "3"], " ", 
                    SuperscriptBox["\[Beta]", "2"]}], "+", 
                   RowBox[{
                   "Y", " ", "z", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"2", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}], 
                    " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "+", 
                   RowBox[{
                    SuperscriptBox["z", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Alpha]", "-", "\[Beta]"}], ")"}], " ", 
                    SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}]}], 
                  ")"}]}], "+", 
                RowBox[{"q", " ", "z", " ", "\[Alpha]", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    SuperscriptBox["Y", "2"], " ", 
                    SuperscriptBox["\[Alpha]", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"6", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                    " ", 
                    SuperscriptBox["\[Beta]", "2"]}], "+", 
                   RowBox[{
                   "Y", " ", "z", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"8", " ", "\[Alpha]"}], "+", 
                    RowBox[{"3", " ", "\[Beta]"}]}], ")"}], " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "+", 
                   RowBox[{"2", " ", 
                    SuperscriptBox["z", "2"], " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Alpha]", "-", "\[Beta]"}], ")"}], " ", 
                    SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}]}], 
                  ")"}]}]}], ")"}]}]}], ")"}]}], ")"}]}]}]}], ")"}], "/", 
     RowBox[{"(", 
      RowBox[{"q", " ", "z", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"64", " ", 
          SuperscriptBox["q", "3"], " ", 
          SuperscriptBox["z", "3"]}], "-", 
         RowBox[{"96", " ", 
          SuperscriptBox["q", "2"], " ", 
          SuperscriptBox["z", "2"], " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
         RowBox[{"48", " ", "q", " ", "z", " ", 
          SuperscriptBox["\[Alpha]", "2"], " ", 
          SuperscriptBox["\[Beta]", "2"]}], "-", 
         RowBox[{
          SuperscriptBox["\[Alpha]", "2"], " ", 
          SuperscriptBox["\[Beta]", "3"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"8", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}]}]}], 
        ")"}]}], ")"}]}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"f\[Sigma]xup", "[", 
    RowBox[{
    "\[Sigma]y_", ",", "\[Alpha]_", ",", "Y_", ",", "z_", ",", "q_", ",", 
     "\[Beta]_", ",", "L_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"64", " ", 
       SuperscriptBox["q", "4"], " ", 
       SuperscriptBox["z", "4"], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"L", " ", "z"}], "+", 
         RowBox[{"Y", " ", "\[Alpha]"}]}], ")"}]}], "+", 
      RowBox[{"8", " ", 
       SuperscriptBox["q", "2"], " ", 
       SuperscriptBox["z", "2"], " ", 
       SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"4", "+", 
            RowBox[{"11", " ", "L"}]}], ")"}], " ", "z", " ", "\[Beta]"}], 
         "+", 
         RowBox[{"7", " ", "Y", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
         RowBox[{"4", " ", "z", " ", "\[Sigma]y"}]}], ")"}]}], "+", 
      RowBox[{
       SuperscriptBox["\[Alpha]", "3"], " ", 
       SuperscriptBox["\[Beta]", "3"], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"\[Beta]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", "\[Alpha]", " ", 
             RowBox[{"(", 
              RowBox[{"z", "+", 
               RowBox[{"2", " ", "L", " ", "z"}], "+", 
               RowBox[{"Y", " ", "\[Alpha]"}]}], ")"}]}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"1", "+", "L"}], ")"}], " ", "z", " ", "\[Beta]"}]}], 
           ")"}]}], "+", 
         RowBox[{"z", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}], " ", 
          "\[Sigma]y"}]}], ")"}]}], "-", 
      RowBox[{"32", " ", 
       SuperscriptBox["q", "3"], " ", 
       SuperscriptBox["z", "3"], " ", "\[Alpha]", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"3", " ", "Y", " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
         RowBox[{"z", " ", 
          RowBox[{"(", 
           RowBox[{"\[Beta]", "+", 
            RowBox[{"4", " ", "L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], 
           ")"}]}]}], ")"}]}], "-", 
      RowBox[{"q", " ", "z", " ", 
       SuperscriptBox["\[Alpha]", "2"], " ", 
       SuperscriptBox["\[Beta]", "2"], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Y", " ", "\[Alpha]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"16", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], " ", 
          "\[Beta]"}], "+", 
         RowBox[{"2", " ", "z", " ", "\[Beta]", " ", 
          RowBox[{"(", 
           RowBox[{"\[Beta]", "+", 
            RowBox[{"L", " ", "\[Beta]"}], "+", "\[Sigma]y"}], ")"}]}], "+", 
         RowBox[{"4", " ", "z", " ", "\[Alpha]", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"3", "+", 
               RowBox[{"7", " ", "L"}]}], ")"}], " ", "\[Beta]"}], "+", 
            RowBox[{"3", " ", "\[Sigma]y"}]}], ")"}]}]}], ")"}]}], "-", 
      RowBox[{"2", " ", "q", " ", "z", " ", 
       RowBox[{"\[Sqrt]", 
        RowBox[{"(", 
         RowBox[{
          FractionBox["1", 
           RowBox[{
            SuperscriptBox["q", "2"], " ", 
            SuperscriptBox["z", "2"]}]], 
          SuperscriptBox["\[Alpha]", "2"], " ", 
          SuperscriptBox["\[Beta]", "2"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", "q", " ", "z"}], "-", 
            RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", 
             SuperscriptBox["L", "2"], " ", 
             SuperscriptBox["z", "2"], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"64", " ", 
                SuperscriptBox["q", "5"], " ", 
                SuperscriptBox["z", "5"]}], "+", 
               RowBox[{"4", " ", 
                SuperscriptBox["q", "3"], " ", 
                SuperscriptBox["z", "3"], " ", "\[Alpha]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"53", " ", "\[Alpha]"}], "-", 
                  RowBox[{"2", " ", "\[Beta]"}]}], ")"}], " ", 
                SuperscriptBox["\[Beta]", "2"]}], "+", 
               RowBox[{"q", " ", "z", " ", 
                SuperscriptBox["\[Alpha]", "3"], " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"24", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                " ", 
                SuperscriptBox["\[Beta]", "4"]}], "-", 
               RowBox[{"2", " ", 
                SuperscriptBox["\[Alpha]", "5"], " ", 
                SuperscriptBox["\[Beta]", "5"]}], "+", 
               RowBox[{"4", " ", 
                SuperscriptBox["q", "4"], " ", 
                SuperscriptBox["z", "4"], " ", "\[Beta]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"-", "48"}], " ", "\[Alpha]"}], "+", "\[Beta]"}], 
                 ")"}]}], "+", 
               RowBox[{"2", " ", 
                SuperscriptBox["q", "2"], " ", 
                SuperscriptBox["z", "2"], " ", 
                SuperscriptBox["\[Alpha]", "2"], " ", 
                SuperscriptBox["\[Beta]", "3"], " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"-", "53"}], " ", "\[Alpha]"}], "+", 
                  RowBox[{"3", " ", "\[Beta]"}]}], ")"}]}]}], ")"}]}], "+", 
            RowBox[{"L", " ", "z", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"4", " ", "q", " ", "z"}], "-", 
               RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"64", " ", 
                SuperscriptBox["q", "4"], " ", "Y", " ", 
                SuperscriptBox["z", "4"], " ", "\[Alpha]"}], "+", 
               RowBox[{"2", " ", 
                SuperscriptBox["q", "2"], " ", 
                SuperscriptBox["z", "2"], " ", "\[Alpha]", " ", "\[Beta]", 
                " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"Y", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"56", " ", "\[Alpha]"}], "+", 
                    RowBox[{"5", " ", "\[Beta]"}]}], ")"}]}], "+", 
                  RowBox[{"2", " ", "z", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"14", " ", "\[Alpha]"}], "-", 
                    RowBox[{"3", " ", "\[Beta]"}]}], ")"}], " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], ")"}]}],
                "-", 
               RowBox[{"2", " ", "q", " ", "z", " ", 
                SuperscriptBox["\[Alpha]", "2"], " ", 
                SuperscriptBox["\[Beta]", "2"], " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"3", " ", "Y", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"6", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}]}],
                   "+", 
                  RowBox[{"2", " ", "z", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"7", " ", "\[Alpha]"}], "-", 
                    RowBox[{"2", " ", "\[Beta]"}]}], ")"}], " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], ")"}]}],
                "+", 
               RowBox[{
                SuperscriptBox["\[Alpha]", "3"], " ", 
                SuperscriptBox["\[Beta]", "3"], " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"Y", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"4", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}]}],
                   "+", 
                  RowBox[{"z", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"4", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                   " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], ")"}]}],
                "-", 
               RowBox[{"8", " ", 
                SuperscriptBox["q", "3"], " ", 
                SuperscriptBox["z", "3"], " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"Y", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"18", " ", "\[Alpha]"}], "+", "\[Beta]"}], 
                    ")"}]}], "+", 
                  RowBox[{"z", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"4", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                   " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}]}], 
                 ")"}]}]}], ")"}]}], "+", 
            RowBox[{
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                 RowBox[{"-", "4"}], " ", "q", " ", "z"}], "+", 
                RowBox[{"\[Alpha]", " ", "\[Beta]"}]}], ")"}], "2"], " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"8", " ", 
                SuperscriptBox["q", "3"], " ", 
                SuperscriptBox["Y", "2"], " ", 
                SuperscriptBox["z", "3"], " ", 
                SuperscriptBox["\[Alpha]", "2"]}], "+", 
               RowBox[{"2", " ", 
                SuperscriptBox["q", "2"], " ", 
                SuperscriptBox["z", "2"], " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   SuperscriptBox["Y", "2"], " ", 
                   SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"-", "6"}], " ", "\[Alpha]"}], "+", "\[Beta]"}], 
                    ")"}]}], "-", 
                  RowBox[{"2", " ", "Y", " ", "z", " ", "\[Alpha]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"2", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}], 
                   " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "+", 
                  RowBox[{
                   SuperscriptBox["z", "2"], " ", 
                   SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}]}], 
                 ")"}]}], "-", 
               RowBox[{
                SuperscriptBox["\[Alpha]", "2"], " ", "\[Beta]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   SuperscriptBox["Y", "2"], " ", 
                   SuperscriptBox["\[Alpha]", "3"], " ", 
                   SuperscriptBox["\[Beta]", "2"]}], "+", 
                  RowBox[{"Y", " ", "z", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"2", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}], 
                   " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "+", 
                  RowBox[{
                   SuperscriptBox["z", "2"], " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Alpha]", "-", "\[Beta]"}], ")"}], " ", 
                   SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}]}], 
                 ")"}]}], "+", 
               RowBox[{"q", " ", "z", " ", "\[Alpha]", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   SuperscriptBox["Y", "2"], " ", 
                   SuperscriptBox["\[Alpha]", "2"], " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"6", " ", "\[Alpha]"}], "-", "\[Beta]"}], ")"}], 
                   " ", 
                   SuperscriptBox["\[Beta]", "2"]}], "+", 
                  RowBox[{"Y", " ", "z", " ", "\[Alpha]", " ", "\[Beta]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"8", " ", "\[Alpha]"}], "+", 
                    RowBox[{"3", " ", "\[Beta]"}]}], ")"}], " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}]}], "+", 
                  RowBox[{"2", " ", 
                   SuperscriptBox["z", "2"], " ", 
                   RowBox[{"(", 
                    RowBox[{"\[Alpha]", "-", "\[Beta]"}], ")"}], " ", 
                   SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"\[Beta]", "+", "\[Sigma]y"}], ")"}], "2"]}]}], 
                 ")"}]}]}], ")"}]}]}], ")"}]}], ")"}]}]}]}], ")"}], "/", 
    RowBox[{"(", 
     RowBox[{"q", " ", "z", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"64", " ", 
         SuperscriptBox["q", "3"], " ", 
         SuperscriptBox["z", "3"]}], "-", 
        RowBox[{"96", " ", 
         SuperscriptBox["q", "2"], " ", 
         SuperscriptBox["z", "2"], " ", "\[Alpha]", " ", "\[Beta]"}], "+", 
        RowBox[{"48", " ", "q", " ", "z", " ", 
         SuperscriptBox["\[Alpha]", "2"], " ", 
         SuperscriptBox["\[Beta]", "2"]}], "-", 
        RowBox[{
         SuperscriptBox["\[Alpha]", "2"], " ", 
         SuperscriptBox["\[Beta]", "3"], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"8", " ", "\[Alpha]"}], "+", "\[Beta]"}], ")"}]}]}], 
       ")"}]}], ")"}]}]}]}]], "Input",
 CellChangeTimes->{{3.962440756376253*^9, 3.962440843596546*^9}},
 CellLabel->"In[52]:=",ExpressionUUID->"3b24fea8-7f11-e844-bfb3-222463853a3f"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"f\[Sigma]xup", "[", 
   RowBox[{
   "\[Sigma]y1", ",", "\[Alpha]1", ",", "Y1", ",", "z1", ",", "q1", ",", 
    "\[Beta]1", ",", "LxA1"}], "]"}], 
  RowBox[{"(*", 
   RowBox[{"\[Sigma]x", " ", "Upper", " ", "Bound"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.9624408581591625`*^9, 3.9624408668934*^9}, {
  3.962529379229025*^9, 3.96252939679994*^9}, {3.962598094672283*^9, 
  3.962598106008882*^9}},
 CellLabel->"In[53]:=",ExpressionUUID->"f8ffa628-3890-844a-83e1-c86943662823"],

Cell[BoxData["0.4880316950224996`"], "Output",
 CellChangeTimes->{
  3.9624408691225224`*^9, {3.962441125225309*^9, 3.9624411434510727`*^9}, 
   3.962441185150156*^9, 3.962441261428835*^9, 3.962441380348423*^9, 
   3.962441422598079*^9, 3.962441538300314*^9, 3.962441614048309*^9, 
   3.962441645017473*^9, 3.962441698765365*^9, 3.9624418036007*^9, 
   3.9624418651441193`*^9, {3.962441919499975*^9, 3.962441945792869*^9}, 
   3.962442194145769*^9, 3.962442417119272*^9, 3.962442713909407*^9, 
   3.9624428809874935`*^9, 3.96244308435864*^9, 3.9624433914275727`*^9, 
   3.96244349703606*^9, 3.962443591203455*^9, 3.962443680789535*^9, 
   3.962443716914226*^9, 3.9625282458311386`*^9, 3.962528322482025*^9, 
   3.962528355315586*^9, 3.962528389076048*^9, 3.962528419842539*^9, 
   3.962528451659662*^9, 3.962528489063364*^9, 3.9625285403669815`*^9, 
   3.962528962565183*^9, 3.962528996112923*^9, 3.962529420893345*^9, 
   3.962530533576828*^9, 3.962530596772394*^9, 3.962531102511345*^9, 
   3.962531162309412*^9, 3.962531407217045*^9, 3.962531439061598*^9, {
   3.962531491175791*^9, 3.962531511271303*^9}, {3.962531544486397*^9, 
   3.962531575201685*^9}, 3.962531633852724*^9, {3.962531703610041*^9, 
   3.962531722152996*^9}, 3.962531779638817*^9, 3.962532109813317*^9, {
   3.9625321419047394`*^9, 3.962532165732256*^9}, 3.962532217997235*^9, 
   3.962532252376436*^9, 3.962532290272984*^9, {3.962532355730289*^9, 
   3.962532376601174*^9}, 3.962532420374563*^9, {3.962532635736194*^9, 
   3.9625326457772655`*^9}, 3.9625326959622765`*^9, 3.9625327383940315`*^9, 
   3.962532851635021*^9, 3.962532944113821*^9, 3.962533019716772*^9, 
   3.962598115952219*^9, 3.96259829785816*^9, 3.9625995893718643`*^9, 
   3.962600088681698*^9, 3.971341518999338*^9},
 CellLabel->"Out[53]=",ExpressionUUID->"4554b00a-c0b9-7c44-b6a0-5a58ab8aaa4a"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1086.857142857143, 539.1428571428571},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"14.0 for Microsoft Windows (64-bit) (2023\:5e7412\:670812\
\:65e5)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"5f8fd964-c762-41ee-a932-cb875a9b5307"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 227, 4, 69, "Chapter",ExpressionUUID->"2f76dbf6-4856-fb4a-b09f-e7ae0b559ca0"],
Cell[CellGroupData[{
Cell[832, 30, 164, 3, 67, "Section",ExpressionUUID->"075a57bb-5672-684f-b1a0-fcf1d8427ebf"],
Cell[999, 35, 879, 21, 48, "Input",ExpressionUUID->"670565b0-2745-7940-9f25-1ce8412da6af"],
Cell[CellGroupData[{
Cell[1903, 60, 762, 15, 28, "Input",ExpressionUUID->"3109ff47-9610-f24b-b339-93226a554405"],
Cell[2668, 77, 756, 16, 32, "Output",ExpressionUUID->"fd18633c-2e83-084a-9f4d-a74f97e793b4"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3461, 98, 928, 18, 28, "Input",ExpressionUUID->"ecf37400-b462-1648-83b7-dc98d91d3794"],
Cell[4392, 118, 836, 17, 32, "Output",ExpressionUUID->"152e8bff-350e-bb4b-ab1c-2c3dfae3c5ef"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5265, 140, 908, 17, 28, "Input",ExpressionUUID->"e73ed0d8-3fe4-ca40-82c3-d5ac7f8f9886"],
Cell[6176, 159, 758, 16, 32, "Output",ExpressionUUID->"6ec3e942-74f5-db4a-99f1-8f8c9cecb68b"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6971, 180, 703, 15, 28, "Input",ExpressionUUID->"99bdc6d7-0dc0-5946-bfdd-8369c67bc491"],
Cell[7677, 197, 753, 17, 32, "Output",ExpressionUUID->"1a5004df-7251-8a4e-b222-4cf13462a382"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8467, 219, 717, 17, 28, "Input",ExpressionUUID->"5d32904f-bbb5-8b4f-bb36-c2fa067cb57b"],
Cell[9187, 238, 699, 14, 32, "Output",ExpressionUUID->"4c00b93b-502b-9f4b-a4ec-24319eb8e863"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9923, 257, 668, 16, 28, "Input",ExpressionUUID->"92b8c5ca-2962-cc49-b1d3-b66eb82809d2"],
Cell[10594, 275, 729, 15, 32, "Output",ExpressionUUID->"314ef302-c208-0942-96a4-07c67cd0bb57"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[11372, 296, 170, 3, 67, "Section",ExpressionUUID->"3aa77a04-0892-7645-acba-1d32f23a05b0"],
Cell[CellGroupData[{
Cell[11567, 303, 1754, 35, 62, "Input",ExpressionUUID->"4aeeeca2-80de-ef43-b694-c81a806f132d"],
Cell[13324, 340, 785, 14, 48, "Output",ExpressionUUID->"852cecc7-adc5-4340-83d6-2e81ee5c8ad2"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14146, 359, 1274, 28, 62, "Input",ExpressionUUID->"d716037c-e532-c942-844f-e2904edaa450"],
Cell[15423, 389, 712, 12, 48, "Output",ExpressionUUID->"5a157eb2-1602-d34d-a08c-ad7534086220"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16172, 406, 1435, 31, 62, "Input",ExpressionUUID->"0f77b9f7-6d6a-d24f-9e59-8212aea8fc7b"],
Cell[17610, 439, 715, 13, 51, "Output",ExpressionUUID->"e704a5ef-27b6-944c-afce-32ebcbbde64e"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18362, 457, 1243, 28, 62, "Input",ExpressionUUID->"fac2a9da-1722-e446-b9f4-1085badaa6db"],
Cell[19608, 487, 618, 11, 51, "Output",ExpressionUUID->"081f3fdd-ed1f-a646-a89b-4078e572e193"]
}, Open  ]],
Cell[CellGroupData[{
Cell[20263, 503, 4040, 87, 233, "Input",ExpressionUUID->"1da6c989-f96d-ae4b-a88d-e89c92f3d6b2"],
Cell[24306, 592, 3169, 80, 89, "Output",ExpressionUUID->"5d8a33c3-7c88-a247-bdf8-2ef28d1f7750"]
}, Open  ]],
Cell[27490, 675, 2564, 52, 162, "Input",ExpressionUUID->"8272f2e7-4d98-de46-9543-3c85c50de20a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[30091, 732, 154, 3, 67, "Section",ExpressionUUID->"58113e13-913c-294c-af9e-4f344f889b47"],
Cell[CellGroupData[{
Cell[30270, 739, 2866, 62, 156, "Input",ExpressionUUID->"287828de-3953-2148-ae4e-5eb78dbcdafa"],
Cell[33139, 803, 1938, 48, 53, "Output",ExpressionUUID->"eec37dfd-6fc2-c942-9f3b-2d3b474828a2"],
Cell[35080, 853, 1958, 52, 50, "Output",ExpressionUUID->"122abe37-9b85-8949-a38d-2c7d7a00de03"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[37087, 911, 173, 3, 67, "Section",ExpressionUUID->"e55b14d4-82f5-b349-828b-ce20e0dafd22"],
Cell[CellGroupData[{
Cell[37285, 918, 2572, 56, 192, "Input",ExpressionUUID->"5e8bcf18-5645-9d45-a88a-b5798abcbb00"],
Cell[39860, 976, 1712, 41, 50, "Output",ExpressionUUID->"df0cdac6-ae92-604a-ba64-7da497281222"]
}, Open  ]],
Cell[CellGroupData[{
Cell[41609, 1022, 1943, 46, 171, "Input",ExpressionUUID->"ffd973ad-bde9-6242-ab8c-b055aca7147c"],
Cell[43555, 1070, 1678, 41, 50, "Output",ExpressionUUID->"3c59924c-fd64-9a41-82e6-b939844d7bf0"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[45282, 1117, 164, 3, 67, "Section",ExpressionUUID->"15f2b005-cb3f-2241-a947-9897c46b1b9a"],
Cell[CellGroupData[{
Cell[45471, 1124, 2755, 62, 159, "Input",ExpressionUUID->"5aacb4e4-5e09-b64d-b3f2-c7f1e9d05f33"],
Cell[48229, 1188, 5074, 133, 175, "Output",ExpressionUUID->"b3f21212-24bd-984a-9afd-849e10c687ac"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[53352, 1327, 174, 3, 67, "Section",ExpressionUUID->"eaf531c3-bdfe-1849-a57a-f69925b30ae8"],
Cell[CellGroupData[{
Cell[53551, 1334, 2527, 57, 140, "Input",ExpressionUUID->"426004d3-426f-744b-8e72-dabc00fae52d"],
Cell[56081, 1393, 2956, 83, 90, "Output",ExpressionUUID->"c35f5ddc-1bae-7240-a7b9-eb4ab2a0d4aa"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[59086, 1482, 274, 4, 67, "Section",ExpressionUUID->"19fb375e-efa7-834e-97ec-10cfe7b7a2ba"],
Cell[CellGroupData[{
Cell[59385, 1490, 463, 8, 28, "Input",ExpressionUUID->"5815f8ef-544d-264f-b39e-3bb52976e630"],
Cell[59851, 1500, 378, 5, 32, "Output",ExpressionUUID->"12850f40-705a-fc42-952d-fd60d00b9927"]
}, Open  ]],
Cell[CellGroupData[{
Cell[60266, 1510, 415, 7, 28, "Input",ExpressionUUID->"7a63be0e-5932-7041-8d96-665aefd17e2a"],
Cell[60684, 1519, 359, 5, 32, "Output",ExpressionUUID->"04e313ed-031f-ed46-ad6e-5bc79eeb9947"]
}, Open  ]],
Cell[CellGroupData[{
Cell[61080, 1529, 515, 9, 28, "Input",ExpressionUUID->"37f57791-5eed-014d-8c2d-5b961029eeae"],
Cell[61598, 1540, 354, 5, 32, "Output",ExpressionUUID->"a1666f0a-40da-cc41-908d-acdb756f1f9b"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[62013, 1552, 173, 3, 69, "Chapter",ExpressionUUID->"31228355-b033-5941-8aa4-15146b7ed2f5"],
Cell[CellGroupData[{
Cell[62211, 1559, 276, 4, 67, "Section",ExpressionUUID->"be53c7ca-f7a6-394e-8ddd-26bff516b8ca"],
Cell[CellGroupData[{
Cell[62512, 1567, 1899, 45, 102, "Input",ExpressionUUID->"f0cbbfc3-6237-8246-b486-bd41e4e20838"],
Cell[64414, 1614, 1924, 44, 50, "Output",ExpressionUUID->"53d859d5-bfcc-ad44-88cf-62f8d5753584"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[66387, 1664, 269, 4, 67, "Section",ExpressionUUID->"cd991f6d-b2a6-5646-b66c-1efca7415d69"],
Cell[CellGroupData[{
Cell[66681, 1672, 5442, 122, 309, "Input",ExpressionUUID->"64219419-dd83-4e49-862c-54ab2a83ce12"],
Cell[72126, 1796, 4285, 106, 173, "Output",ExpressionUUID->"a551a93a-f316-c548-b74d-cf357c984c9e"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[76460, 1908, 255, 4, 67, "Section",ExpressionUUID->"195018d1-e59b-624f-996a-cebc149953f4"],
Cell[76718, 1914, 572, 17, 28, "Input",ExpressionUUID->"d13025f6-611f-8a40-a1bd-90d159f1cb41"],
Cell[CellGroupData[{
Cell[77315, 1935, 1286, 31, 83, "Input",ExpressionUUID->"ca6399b8-d9f7-cc43-9c6c-5580c076bcc3"],
Cell[78604, 1968, 3069, 71, 94, "Output",ExpressionUUID->"c1b0a2f9-a635-6d47-9409-07fa31e3a155"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[81722, 2045, 155, 3, 67, "Section",ExpressionUUID->"e8996b59-8377-1745-b19d-5cbad812efca"],
Cell[CellGroupData[{
Cell[81902, 2052, 1186, 30, 62, "Input",ExpressionUUID->"56c43505-5841-714b-a949-3ba28d803d6a"],
Cell[83091, 2084, 1723, 39, 50, "Output",ExpressionUUID->"3de31723-4264-7748-88bf-f9b189369f47"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[84863, 2129, 345, 5, 67, "Section",ExpressionUUID->"09599678-18b0-0446-ae13-db4e325c73bc"],
Cell[CellGroupData[{
Cell[85233, 2138, 2276, 51, 156, "Input",ExpressionUUID->"31f1adab-b60d-ac44-9856-6da97ad627db"],
Cell[87512, 2191, 16367, 392, 317, "Output",ExpressionUUID->"90cc936a-c9b5-6344-b64e-68c625c9e60d"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[103928, 2589, 235, 4, 67, "Section",ExpressionUUID->"affa3178-3bb7-7b44-b932-16da0c4dbf2b"],
Cell[CellGroupData[{
Cell[104188, 2597, 804, 21, 43, "Input",ExpressionUUID->"ec5783c5-fa28-334f-83f0-c5e0f7ccfcbf"],
Cell[104995, 2620, 417, 10, 28, "Message",ExpressionUUID->"567307af-5a62-7b49-9e73-9b84b3b850f5"],
Cell[105415, 2632, 13613, 332, 234, "Output",ExpressionUUID->"6bdae327-aa66-2c49-83cf-a9b64c6e5cfa"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[119089, 2971, 173, 3, 69, "Chapter",ExpressionUUID->"1a2d8b0b-902e-e44e-ab23-ff4488f72778"],
Cell[CellGroupData[{
Cell[119287, 2978, 403, 8, 67, "Section",ExpressionUUID->"304c1d60-8da9-c14b-b4f8-da5fe23ed6ba"],
Cell[CellGroupData[{
Cell[119715, 2990, 6417, 138, 271, "Input",ExpressionUUID->"3ace2da8-68f0-8d42-a9ae-1f9137c2cb87"],
Cell[126135, 3130, 4691, 137, 96, "Output",ExpressionUUID->"cf9dd6cf-8edd-b141-8596-780e65f5b6a9"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[130875, 3273, 229, 4, 67, "Section",ExpressionUUID->"0b4a5684-4873-3d44-bda0-f14e89ac48a7"],
Cell[CellGroupData[{
Cell[131129, 3281, 4787, 116, 194, "Input",ExpressionUUID->"2ebbe932-4037-4f44-ba81-65c14d504b23"],
Cell[135919, 3399, 1101, 29, 50, "Output",ExpressionUUID->"1cbb047e-f000-684d-93f2-01bfecb71da4"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[137069, 3434, 247, 4, 67, "Section",ExpressionUUID->"40806356-fd71-5b41-a5d6-56bd84606337"],
Cell[CellGroupData[{
Cell[137341, 3442, 3789, 100, 105, "Input",ExpressionUUID->"747f4495-10dc-6349-810f-acffd2fa190d"],
Cell[141133, 3544, 3206, 89, 93, "Output",ExpressionUUID->"345ca96e-be52-af47-9452-f980d3d1d221"],
Cell[144342, 3635, 3426, 92, 99, "Output",ExpressionUUID->"7e48b32b-1706-2347-816e-b144472b921b"],
Cell[147771, 3729, 3318, 97, 90, "Output",ExpressionUUID->"3a0560b4-1d04-8e4d-9555-fa7cb7aca224"],
Cell[151092, 3828, 3555, 102, 99, "Output",ExpressionUUID->"f2bc6e6b-c162-f446-88f2-297252e4cf48"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[154708, 3937, 161, 3, 69, "Chapter",ExpressionUUID->"4945309f-a7de-5f42-9796-cd1dd11225cd"],
Cell[CellGroupData[{
Cell[154894, 3944, 228, 4, 67, "Section",ExpressionUUID->"abf52c42-3e69-8b41-bb64-f9463486e6d8"],
Cell[155125, 3950, 3310, 51, 28, "Input",ExpressionUUID->"5fb78389-c9f4-6c45-914d-179ce1329433"],
Cell[CellGroupData[{
Cell[158460, 4005, 1917, 47, 121, "Input",ExpressionUUID->"90ff39bb-5d09-c644-b727-ea98751d0791"],
Cell[160380, 4054, 1329, 20, 32, "Output",ExpressionUUID->"6e9bd15f-2e4b-8841-9c12-d07ed6d5c49a"]
}, Open  ]],
Cell[161724, 4077, 2200, 64, 85, "Input",ExpressionUUID->"002a7799-b916-2249-ab40-b7777a612d04"],
Cell[CellGroupData[{
Cell[163949, 4145, 497, 11, 28, "Input",ExpressionUUID->"9141a280-bbf2-7040-8c29-1a18e9f8c7bb"],
Cell[164449, 4158, 1881, 26, 32, "Output",ExpressionUUID->"c02a5c13-1f18-5346-820e-8dedf1e79bf1"]
}, Open  ]],
Cell[CellGroupData[{
Cell[166367, 4189, 499, 11, 28, "Input",ExpressionUUID->"244b78b3-f103-da4a-85ff-469d432fcbf2"],
Cell[166869, 4202, 1863, 26, 32, "Output",ExpressionUUID->"de9d5bc4-08fe-1344-b03c-6161d2a9d247"]
}, Open  ]],
Cell[168747, 4231, 482, 11, 48, "Input",ExpressionUUID->"d8f4e21a-0f1e-1648-8a0a-e09550d07089"],
Cell[CellGroupData[{
Cell[169254, 4246, 420, 9, 28, "Input",ExpressionUUID->"90caf400-2533-b940-9809-a224eeef115b"],
Cell[169677, 4257, 1700, 24, 32, "Output",ExpressionUUID->"f8d7c523-ac41-9f4e-a6a2-eefc5d4e81a7"]
}, Open  ]],
Cell[171392, 4284, 28368, 673, 644, "Input",ExpressionUUID->"3b24fea8-7f11-e844-bfb3-222463853a3f"],
Cell[CellGroupData[{
Cell[199785, 4961, 513, 11, 28, "Input",ExpressionUUID->"f8ffa628-3890-844a-83e1-c86943662823"],
Cell[200301, 4974, 1837, 26, 54, "Output",ExpressionUUID->"4554b00a-c0b9-7c44-b6a0-5a58ab8aaa4a"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)

